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Computer Modeling of Flying Star Ballistics 


Dayu Ding, Morimasa Higaki, Yozo Ooki and Tadao Yoshida* 


Ashikaga Institute of Technology, 
268-1 Omae-cho, Ashikaga-shi, Tochigi 326-8558, KAPAN 
Tel: +81-284-62-0605 
Fax: +81-284-62-0976 
E-mail yoshida@ashitech.ac.jp 
*To whom all correspondence should be addressed. 


Abstract: The burning time of stationary and flying fireworks stars was measured, the trajectories of flying 
stars were observed and the results were analyzed in this work. It was found that the difference in burning 
time between the stationary and flying burning stars was dependent on the kind of star. Modeling of flying 
star ballistics was applied to the trajectory of stars with shorter burning times and was found valid for this 


case. 


Keywords. fireworks, burning time, exterior ballistics, modeling 


Introduction 


Firework stars with short burning times are used 
for warimono aerial shells giving spherical fire 
flowers, katamono aerial shells giving various burst 
shapes, Roman candles, and so on. The ballistics 
of stars with short burning times are important for 
designing warimono and katamono shells, Roman 
candles and others. 


The present authors have suggested a computer 
model for estimating the trajectory of the burning 
star expelled from a mortar using a star with longer 
burning time.' In the present work, the model is 
applied to a flying star with a shorter burning time 
and the results are analyzed. 


Computer modeling of aerial shells has been done 
by Kosanke and Kosanke,” and Mercer.’ In the case 
of an aerial shell, the drag coefficient (Cp) can be 
assumed constant before the bursting of the shell 
in the air. But, in the case of a burning star, the 
above assumption is not valid, and an approximate 
approach was adopted in our previous work.! 


Experimental 
Materials 


The silver peony stars for no. 2, 2.5, 3, 4 and 5 
shells (a Japanese no. 2.5 shell corresponds to a 
Western 3 inch shell), lifting charge and electric 
matches were supplied by Sunaga Fireworks Co. 
Ltd. at Ashikaga City. 
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Grain black powder made by Nippon Kayaku 
Co. Ltd. was used as the lifting charge in our 
experiments. The standard and particle distribution 
of the lifting charge were described in our previous 
paper.” 


Apparatus 


The three kinds of mortar used for firing stars 
were made of steel, and the inner diameters were 
12mm, 15mm and 20mm, and depths were 
216 mm, 270 mm and 360 mm, respectively. 


The burning time and trajectories of stars were 
measured and recorded using a high-speed video 
camera (Phantom VR-V4.2). 


Procedure 


In stationary burning experiments, a star was 
placed on a heat resistant board and ignited by a 
torch. The stationary burning time of the star was 
determined using a high-speed video camera. 


In star shooting experiments, the mortar was set 
on the ground vertically. An electric match was 
placed on the bottom of the mortar, the lifting 
charge was poured into the mortar from the 
muzzle, and then a star was placed on the lifting 
charge. The electric match was ignited by turning 
on an electric current. The star was shot into the 
air after the lifting charge burned. 


The trajectory of the star was recorded by the 
high-speed camera. Each frame of the video was 
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reproduced on a video screen, and the burning 
time and the relationship between flying time and 
flying height of the star were determined. 


Computer modeling of flying star 
ballistics’ 


Theory 


The motion of a burning star expelled vertically in 
the air can be expressed as follows. 


du = g 3D air : Cy : | U 

dt ADs D star ( 1) 
or 

du =-g 3P air .K |u| U (2) 
dt 4D star 


Here, w, t, 25 Pair Pstars star and Cp are star velocity, 
time, acceleration of gravity, air density, star 
density, diameter of the star and drag coefficient of 
air, respectively. The velocity u is positive when 
the star moves upward. 


And 


K= (3) 


star 


It is assumed that the mass and cross sectional area 
of a burning star change but the density of the star 
does not change with time. Therefore, the second 
term of the right hand side of equation (1) for the 
air drag is proportional to the air drag coefficient 
Cp and inversely proportional to the diameter of 
the star Dg, aS expressed in equation (3). 


K is obtained by the step by step calculation of the 
increment of the trajectory using equation (2). The 
diameter D,,, of the star is a function of the flying 
time of the star and is calculated from the linear 
burning rate of the star: 


D star D star0 2r, staré (4) 


Here, Detaro aNd Par are the initial diameter and the 
linear burning rate of the star, respectively. 
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Then, Cp is calculated from K and D,,, using 
equations (3) and (4). It was found that Cp is 
nearly a linear function of time in the earlier stages 
of the trajectory and a scattered complex function 
of time in the later stages.’ It was also found that 
in the low velocity range Cp of firework shells has 
little effect on the air drag.* 


In our modeling of flying star ballistics, Cp is 
approximated as the linear function of time and 
expressed in equation (5): 


Cp =att+b (5) 


aand b in equation (5) are determined from the star 
shooting experiments and are used for estimating 
the trajectory of a burning star. 


Results and Discussion 
Stationary and flying burning times of stars 


The mean burning time of the stationary and 
flying stars is listed in Table 1. In the previous 
experiments, ° t/t, were much larger than 1.0, but 
in this experiment, values of t/t, were nearly 1.0. 
The differences may be attributable to the nature 
of the stars used. 


Experimental and calculated results of the 
trajectories of stars 


The experimental and calculated results of the 
trajectories of stars are listed in Table 2. The 
relative standard deviations of observed data were 
below 10% except for the initial velocity of the 
stars. The absolute values of the relative standard 
deviation of a are large because the mean values of 
a are very small. In this case the relative standard 
deviation has less meaning. 


Fit of the calculated to the observed trajectory 


All of the calculated trajectories were fitted to the 
observed ones. Examples are shown in Figure 1. 


Effect of the kind of stars on a and b 


In our previous work,’ different values of a and 
b in equation (5) were obtained with the silver 
crown stars for no. 4 and no. 5 shells compared to 
those with silver peony stars in this work as listed 
in Table 3. The burning behavior may be different 
with different kinds of stars. 
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Table 1 Mean burning times of the stationary and flying stars. 


Stationary Flying 
Bu Pesaro time(s) sp RSD time (1) SD RSD ts 
(mm) (ms) oe (ms) 
Silver peony for 
No. 2 shell 10.27 2608 163 0.06 10.33 2751 126 0.05 1.05 
No. 2.5 shell 11.19 2833 42 0.01 11.29 3311 197 0.06 1.17 
No. 3 shell 12.13 3340 115 0.03 12.38 3389 318 0.09 1.01 
No. 4 shell 13.30 3659 102 0.03 13.37 3464 153 0.04 0.95 
No. 5 shell 14.91 3809 828 0.22 15.01 3967 138 0.03 1.04 
No. 5 shell* 14.38 3822 195 0.05 1.00 
Silver crown for 
No. 4 shell 16.57 4198 0.04 16.84 5701 0.05 1A 
No. 5 shell 17.64 4646 0.04 17.62 6761 0.03 15° 
Blue peony, silver peony and silver crown for 
No. 2-6 shells 16° 


Halt the ignition promoter was scraped OH and covered by an inhibitor. SD 1s standard deviation, and RSD (=SD/ 


mean) is relative standard deviation. 


Effect of the size and initial velocity of stars on 
a and b 


The calculated results for stars from no. 2 (smallest) 
and no. 5 (largest) shells are shown in Figure 2. 
Values of Cp for the smallest star from a no. 2 shell 
increased with time and Cp) of the largest from a 
no. 5 shell decreased with time. b in equation (5) 
was 0.46—0.53 with the smallest star from a no. 2 
shell, while it was 0.55—0.63 for the largest star 
from a no. 5 shell. 


Plots of a and b against the initial size and velocity 
of the stars are shown in Figure 3. The scatter 
of a of the stars for a no. 2 shell (runs 1-5) is 
abnormally large and the value for run 19 is also 
abnormally small. These data were excluded from 
the statistical consideration. 


The correlations of a and b with Dga.9 and uy do not 
appear significant. Therefore, the mean values of 
aand b calculated from original data were —0.010 
and 0.57, respectively. 


The estimated trajectories of the burning silver 
peony stars using the above values of a and b are 
plotted against time in Figure 4 along with the 
observed trajectories. 
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Table 2 Experimental and calculated results on the trajectories of the silver peony stars. 
aa a CT an a 771 2; 


oa oo a Shell va i an (ms!) ee ig (mm : a 
(mm) (g) (ms) s) 
1 12 0.5 No.2 1.12 10.2 81 2.12 0.069 0.526 
2 12 0.5 No.2 1.05 10.2 61 2772 1.95 0.056 0.455 
3 12 0.5 No.2 0.94 9.8 61 2736 1.86 0.036 0.492 
+ 12 0.5 No.2 1.16 10.9 62 2594 2.12 —0.001 0.528 
> 12 0.5 No.2 1.18 10.5 85 2901 1.87 0.036 0.512 
Mean 1.09 10.3 70 2751 1.99 0.039** 0.503 
SD 0.10 0.4 12 126 0.13 0.026 0.030 
RSD 0.09 0.04 0.17 0.05 0.07 0.669 0.060 
6 12 0.5 No. 2.5 1.46 11.3 87 3041 1.85  —0.001 0.490 
7 12 0.5 No. 2.5 1.45 11.3 92 3478 1.62 —0.003 0.510 
8 12 0.5 No. 2.5 1.40 11.2 100 3436 1.63 —0.006 0.571 
9 12 0.5 No. 2.5 1.52 11.7 80 3288 1.78 -0.012 0.531 
10 12 0.5 No. 2.5 1.36 11.0 86 1.74 —0.004 0.530 
Mean 1.44 11.3 89 3311 1.73. —0.005 0.526 
SD 0.06 0.3 8 197 0.10 0.004 0.030 
RSD 0.04 0.02 0.09 0.06 0.06 —0.835 0.057 
11 15 1 No. 3 2.02 12.6 76 2.30 0.000 0.610 
12 15 1 No. 3 1.95 12.5 87 3808 1.85 —0.009 0.561 
13 15 1 No. 3 1.90 12.5 93 3110 2.10 0.000 0.490 
14 15 1 No. 3 1.89 12.3 95 3462 2.15 0.000 0.530 
15 15 1 No. 3 1.85 12.1 94 3176 2.04 —0.008 0.571 
Mean 1.92 12.4 89 3389 2.09 —0.004 0.552 
SD 0.06 0.2 8 318 0.16 0.005 0.045 
RSD 0.03 0.01 0.09 0.09 0.08 —1.327 0.082 
16 15 1 No. 4 2.43 13.1 111 3300 2.09 —0.033 0.624 
17 15 1 No. 4 2.61 13.6 86 3604 1.99 -0.014 0.612 
18 15 1 No. 4 2.49 13.6 102 2.13 —0.005 0.601 
19 15 1 No. 4 2.38 13.3 108 3488 2.15  —0.071* 0.595 
20 15 1 No. 4 2.34 13.3 103 1.84 —0.003 0.590 
Mean 2.45 13.4 102 3464 2.04 —0.014 0.604 
SD 0.10 0.2 10 153 0.13 0.028 0.014 
RSD 0.04 0.02 0.10 0.04 0.06 —1.110 0.023 
21 20 2 No. 5 3.46 15.2 79 3950 2.20 —0.013 0.612 
22 20 2 No. 5 3.36 14.7 91 3916 2.01 —0.015 0.592 
23 20 2 No. 5 3.58 15.5 115 4198 2.01 —0.014 0.611 
24 20 2 No. 5 3.34 14.7 88 3826 2.03 —0.024 0.635 
25 20 2 No. 5 3.46 15.1 92 3944 2.04 —0.019 0.553 
Mean 3.44 15.0 93 3967 2.06 -—0.017 0.600 
SD 0.10 0.3 13 138 0.08 0.005 0.031 
RSD 0.03 0.02 0.14 0.03 0.04 —0.269 0.051 
Overall mean —0.010 0.560 


* These data are omitted from the mean calculation. ** These data are omitted from the overall 
mean calculation. 
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Table 3 Values of a and b in equation (5) for different stars. 


Mean diameter 


Kind of star 
(mm) 
Silver crown star for 
No. 4 shell 16.8 
No. 5 shell 17.5 
Silver peony star for 
No. 4 shell 13.4 
No. 5 shell 15.0 


Height(m) 


0 0.5 1 1. 


Time(s) 


5 


Observed 
Calculated 


2 


(a) Star for no.2 shell (run1) 


Height(m) 
oo 
ras) 


0 1 
Time(s) 


2 


Observed 
Calculated 


(©) Star for no.3 shell (run 11) 


2.5 


Mean uy 


159 
167 


102 


Height(m) 


Height(m) 


93 


(mm s ') 


+ "Observed 
Calculated 


0 1 2 3 
Time(s) 
(b) Star for no.2.5 shell (run 6) 


Observed 
Calculated 


0 1 2 3 
Time(s) 
(d) Star for no.4 shell (run 16) 


Observed 
—— Calculated 


Time(s) 
(e) Star for no.5 shell (run 21) 


Figure 1 Examples of the observed and calculated trajectories of burning stars. 
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Height(m) 


0 1 2 3 7 1 7 3 0 1 2 3 
Time(s Time(s) Time(s) 
(a) K of stars for no.2 (b) Cp of stars for no.2 (c) Trajectories of stars for no.2 
500 0.7 80 
7 yun 21 0.6 70 
400 7 ¢ run 22 ° - 60 
4 yun 23 all 7 ‘Ee 
“> 300 | x run 24 A 0.4 Can B50 
cS X yun 25 By & 63 5 yun 21 4 40 
200 Be ° [ @ run 22 2 30 
0.2 4 yun 23 
100 fF X run 24 20 
0.1 fF x ym 25 10 
0 1 i ! 0 i 7 ; 
0 Hl 2 3 4 0 1 2 3 4 0 i 2 3 4 
Time(s Time(s) Time(s) 
(d) K of stars for no.5 (e) Cp of stars for no.5 (f) Trajectories of stars for no.5 


Figure 2 Calculated K, Cp and flying height vs. flying time for no. 2 and no. 5 shells. 


0.08 —_—_————————————— 0.08 Lf e 
0.06 | M star for no.2 - 0.06 Ff - @ star for no.2 
® star for no.2.5 ® star for no.2.5 
0.04 | Astar for no.3 LL 0.04 P = Astar for no.3 
0.02 + O star for no.4 0.02 Ff O star for no.4 
& goo L © star for no.5 Ac 30.00 | ma © star for no.5 
-0.02 ta & -0.02 + of - 
-0.04 LE -0.04 i 
-0.06 Ff -0.06 F 
o im] 
-0.08 -0.08 
0 5 10 15 20 0 50 100 150 200 
D staro (mm) uo(ms ) 
(a) Plot of a vs .D staro (b) Plot of a vs. uo 
0.7 7 
0.6 + soy prf'o 
i | 
0.5 7 : 
0.4 7 
— 2 
™ star for no.2 0.3, star for no.2 
® star for no.2.5 02 | ® star for no.2.5 
A star for no.3 : A star for no.3 
O star for no.4 0.1 | Ustar for no.4 
© star for no.5 © star for no.5 
0.0 
0) 5 10 15 20 0 50 ( 4100 150 
D staro (mm) ONE 
(c) Plot of b vs. D staro (d) Plot of b vs. uo 


Figure 3 Plots of a and b vs. the initial size and velocity of burning silver peony stars. 
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xq 20 xX rung x run 14 
X run 10 run 15 
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0 0 it | 
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Time(s) Time(s) 
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(c) Stars for no.4 shell (d) Stars for no.5 shell 


Figure 4 Estimated and observed trajectories of burning stars. 


References 


1 Y. Ooki, D. Ding, M. Higaki and T. Yoshida, 
“Trajectory of Stars with Long Burning 
Time”, Journal of Pyrotechnics, in press. 


2 K.L. Kosanke and B. J. Kosanke, 
“Computer Modeling of Aerial Shell 
Ballistics ”, Pyrotechnica, XIV, 1992, p. 2. 


3. J. E. Mercer, “Thermodynamics of Black 
Powder and Aerodynamics of Propelled 
Aerial Shells”, Journal of Pyrotechnics, 
Issue 16, Winter 2002, p. 37. 


4 — Y. Ooki, D. Ding, M. Higaki and T. Yoshida, 
“Air resistance of Spherical Fireworks 
Shells”, Science and Technology of 
Energetic Materials, in press. 


5D. Ding, M. Higaki, Y. Ooki and T. Yoshida, 
“Pressure in a Mortar and Estimation of 
Muzzle Velocity of Expelled Stars”, Journal 
of Pyrotechnics, Issue 22, Winter 2005, 
p. 50. 


6 Y. Ooki, D. Ding, M. Higaki and T. Yoshida, 
“Burning and Air Resistance of Fireworks 
Stars”, Science and Technology of Energetic 
Materials, in press. 


Journal of Pyrotechnics, Issue 23, Summer 2006 Page 9 


Temperature measurements within the luminous region 
of a burning Ba(NO,)./Al mixture 


P. J. Disimile, R. Prasad and N. Toy 
UC-FEST, Department of Aerospace Engineering, University of Cincinnati, Ohio 45221, USA 


Abstract: Knowledge of the local temperature field associated with a pyrotechnic event has numerous 
implications, especially in the field of safety and survivability. These implications involve the development 
of sensors that are capable of detecting pyrotechnic events and that are used in part to eliminate or reduce 
a fire hazard. However, in order to be able to predict a possible fire scenario from a pyrotechnic event the 
temperature distributions and the thermal heat transfer are prerequisites. This experimental study discusses 
the temperature measurement methodology required to evaluate the transient temperatures associated with 
a small, commercially available, pyrotechnic device. Furthermore, the temperature distribution close to 
the surface of two devices, one commercial, the other fabricated, has been obtained, and shows that the 


temperature distribution away from the event is not uniform. 


Keywords: Zemperature distribution, thermocouples, pyrotechnic facility 


Introduction 


The use of pyrotechnic compositions can be found 
in numerous fields like rocket propellants, highway 
flares, entertainment and other high-energy 
applications. Due to the release of a large amount 
of thermal energy, the burning composition and the 
events occurring around a pyrotechnic device are 
saturated by the intensity of emitted light. Hence 
a uniform region of light is observed which gives 
the impression of a fireball of constant temperature 
throughout. This is even more prominent when 
considering transient pyrotechnic events such as 
in the use of fireworks and incendiary devices. 
The measurement of the high temperatures that 
are produced by such near instantaneous chemical 
reaction is difficult to achieve and challenging 
to understand, and the ability to measure the 
temperature of a dynamic and transitory thermal 
field is far from trivial. 


In order to be able to measure the expected high 
temperatures within the flash/fireball of a short 
duration pyrotechnic event a number of different 
techniques were considered. Of these, the most 
notable that have been used are optical pyrometers 
whereby the brightness of a flame has been comp- 
ared to the brightness of an incandescent filament 
in order to determine the flame temperature.' 
Other types of pyrometers have also been used, but 
have not been successful, for example, cinephoto 
pyrometer, photoelectric photometer and a 
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color photometer.’ In comparison, temperature 
measurements of complex combustion within 
non-uniform temperature zones have been used 
successfully using a line-reversal methodology. 
However, this method only provides intermediate 
temperature distributions within the various zones 
of a flame and is dependent upon the emission 
of characteristic spectral lines from the flame as 
viewed by a spectroscope.” However, even in this 
latter method, the temperature measurements only 
accounted for the average temperature across a 
flame region and the variations within the inner 
zones were unaccounted for. Although temperature 
fields with relation to flame height have been 
measured using a cinephoto pyrometer,' detailed 
temperature distributions were unattainable since 
this method could not provide measurements at 
intervals less than 10 mm. 


In order to satisfy the goals of the present study, 
it was considered necessary to adopt an approach 
of using high-speed thermocouples to capture 
the temperature distribution. Such temperature 
measurements have been performed earlier at 
a single fixed position in a closed system on 
different pyrotechnic materials, notably Sb/ 
KMnO, (antimony/potassium permanganate),’ 
and Pd/Al (palladium/aluminum) mixtures.* 
Temperature profiles have also been measured 
using W/Re (tungsten/rhenitum) thermocouples 
for temperatures above 2000°C in Mo/KCIO, 
(molybdenum/potassium perchlorate) material,° 
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and W/KCIO,/BaCrO,  (tungsten/potassium 
perchlorate/barium chromate) material,° but again 
only at a fixed position in the system. Although 
bare bead W/Re type thermocouples are capable 
of measuring temperatures in excess of 2000 °C, 
they are very prone to oxidation that can lead 
to large errors and have to be used within inert 
environments to prevent the tungsten from 
oxidizing. 

The ability to measure and to analyze the dynamic 
thermal events surrounding a short duration 
pyrotechnic is problematical but it was not 
considered impractical if the temperature range of 
the apparatus could be obtained and if the thermal 
response was fast enough to record the event. Given 
these complications, it was decided to examine the 
capability of miniature high-speed thermocouples 
to observe the temperature distribution close to 
the surface of a pyrotechnic event. Two types 
of thermocouple were initially chosen for this 
study, an R- and a K-type which have the ability 
to measure maximum temperatures up to 1750 °C 
and 1250°C respectively. A pyrotechnic event 
was created using, in the first instance, a common 
sparkler, and secondly a laboratory made device 
containing a chemical mixture similar to that of a 
commercial pyrotechnic material (CPM). 


Experimental Strategy 


Since thermocouples are available in different 
combinations of metals or calibrations it was 
necessary to select ones that were favorable to this 
study. The four most common types are J, K, T and 
E, with each type having a different temperature 


Table 1 Effect of different thermocouple constructions. 


Response 


range and environmental usage, although the 
maximum temperature that each may record 
depends on the diameter of the wire used in the 
thermocouple. Given the harsh environment of the 
present studies, the following criteria were used in 
selecting a suitable thermocouple: 


(1) Temperature range: The operating range 
of a thermocouple is the temperature 
range over which the thermocouple will 
perform satisfactorily, with negligible error 
in the output signal. For the current set of 
experiments, a temperature range with a high 
upper measuring limit was desirable. 


(2) Thermocouple junction selection: Each 
thermocouple mustutilize a measuring junction 
and a reference junction at two different 
temperatures. The measuring junction is 
generally at the higher of the two temperatures 
and the reference junction is at ambient. 
The measuring junction is placed near or on 
whatever is to be measured and the reference 
junction is connected either to a controller or 
a temperature indicator. Different kinds of 
measuring junctions are used with respect to 
measuring requirements as shown in Table 1. 
From this table it may be observed that for 
the present study an exposed bead weld, with 
a low thermal mass and corresponding fast 
response time, would be the most acceptable 
geometry for a thermocouple. 


(3) Response time using first order response 
criteria: The time constant is the time required 
for the measured temperature to reach 63.2% 


Type of junction fies Advantages Disadvantages 
Reliabl : : 
Ungrounded Slow ean a mugged Sluggish response time 
construction 
Sheathed . 
Grounded Normal wees = ascetics) Noise injection 
conductive metallic sheaths 
Bead weld Fast Low thermal mass Prone to damage in a corrosive 
increases response time environment 
Exposed — 
Butt weld Fastest Useful in high speed Corrosive failure and physical 
measurements or mechanical damage 
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of a step change in the temperature of the 
surrounding media. Five time constants are 
required for the sensor to approach 100% of 
the step change value. Typical response time 
for thermocouples range from milliseconds 
to seconds, depending on the size of the 
thermocouple, the fluid thermal conditions 
and the junction employed. 


(4) Operating environment: The operating 
temperature and external environment can 
affect the performance of the thermocouple. 
The present test conditions required 
thermocouples whose output signals were 
not affected by any particulate from the flash/ 
fireball scenario. 


In order to record the temperature measurement 
from the selected thermocouple a dedicated 
acquisition system was used. This consisted of 
a 16-channel analog to digital data acquisition 
system with 12-bitresolution (National Instruments 
PXI - E6040E) that was connected to a stand- 
alone personal computer. LabVIEW software 
controlled the data acquisition program in such a 
way that different types of thermocouples could 
be accommodated using a single input channel. 
This approach was adopted so that the two types 
of thermocouple (R & K) could be tested through 
the same channel. Since the pyrotechnic event has 
to be captured within a short time period the data 
acquisition rate was set at 100 samples per second, 
with more than 1500 samples being recorded, 
depending on the test procedure. 


It should be noted that initially a C-type 
thermocouple was also selected along with the R- 
and K-types because of their ability to withstand 
temperatures up to 2320 °C. The material used 
in the construction of these different types of 
thermocouples consists of two wires of: 


(a) platinum and the other platinum/rhodium 
alloy for the R-type, 


(b) tungsten and rhenium alloy for the C-type, 
and 


(c) chromel (a nickel chromium alloy) and 
alumel (a combination of nickel, aluminum, 
manganese and silicon) for the K-type. 


However, from an initial investigation on the usage 
ofa0.217 mm (36-gauge) R-type thermocouple for 
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measuring the surface temperature of a common 
sparkler both the C- and R-type thermocouples 
were discarded in favour of the K-type. The 
reason for this decision may be considered from 
the following study where a 40mm length of 
a standard 2.8mm diameter sparkler has been 
mounted horizontally on an aluminum support. 
A second aluminum plate supports an R-type 
thermocouple such that the thermocouple bead 
is in contact with the surface of the sparkler, 
Figure 1. 


However, it should be noted that once a sparkler 
is ignited, Figure 2, the visible sparkling is 
created by the release of small iron particles that 
are distributed throughout its body and which 
may also be observed as physical protuberances 


Mounting plate 


> 40-mm of ‘sparkler’ 


R-type thermocouple 


Figure 1 An R-type thermocouple touching the 
surface of a horizontally held sparkler. 


lron particles 


" 
j 


Advancing flame front 


Figure 2 An advancing flame front passes over 
the R-type thermocouple bead. 
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Figure 3 Surface temperature of a sparkler using 
an R-type thermocouple. 


on its surface. In order that the thermocouple 
sensor could be located accurately against the 
sparkler surface these protrusions were removed 
before mounting the specimen in the supporting 
plate. As the flame front passes the thermocouple 
bead, the local surface temperature was recorded 
and a typical temperature-time plot is shown in 
Figure 3. 


Here it is shown that as the flame front approaches 
the bead there is a rapid increase in temperature 
followed by series of temperature perturbations, 
reaching a maximum temperature of approximately 
1575° C, before falling rapidly once the flame front 
has passed. These perturbations were considered 
related to the release of the small iron particles that, 
on their own, do not add to the heat generation but 
are likely to cause a local temperature fluctuation 
on their discharge. After considering the potential 
for large temperature errors to occur due to the 
contamination of the measuring bead at these high 
temperatures where the hot gases produced by 
the reaction cannot be controlled, a decision was 
made to discard their use in favour of the K-type 
thermocouple. 


Similarly, a 0.217 mm (36-gauge) bare bead, 
unsheathed, C-type thermocouple was also 
examined, and not selected. Every attempt to 
sheath the bare wires resulted in damage to the 
thermocouple near the bead, due to the brittle 
nature and associated difficulty in handling of the 
wires. Furthermore, the tungsten wire required an 
inert atmosphere, since it also quickly oxidized at 
high temperatures. 
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The thermocouple eventually chosen was a 
0.0787 mm diameter (40 gauge) K-type with an 
exposed bead weld junction. This thermocouple 
consists of two wires, which are insulated from 
each other except in the region of the bead, 
placed inside a sheath of 0.0762 mm thickness. 
The positive wire is made of chromel, which is a 
composition of 90% nickel and 10% chromium, 
and the negative wire is made of alumel, a 
composition of 95% nickel, 2% aluminum, 2% 
manganese and 1% silicon. The advantages of 
using the K-type thermocouples are: 


(1) They have the highest temperature range 
(-200°C to 1250°C) among the most 
commonly used types of thermocouples. 


(2) They could be fabricated within the laboratory. 
This was achieved by using pre-insulated wire 
and removing a 15 mm length of insulation 
from each. The exposed wires were then 
twisted and welded into a small bead in a 
thermocouple welder. 


When the welded connection (bead) of a 
thermocouple is heated, a voltage across the two 
junctions is produced. A polynomial equation is 
then used to convert the thermocouple voltage 
(E) to a temperature (7/°C) over a wide range of 
temperatures and is given by the equation based 
on the International Temperature Scale (ITS—90) 
standard: 


N . 
(ENO) as 
i=0 


where the coefficients, d; are as given by the 
National Bureau Standards in Table 2 below for 
K-type thermocouples.’ 


Time response evaluations of the K-type 
thermocouple 


In order to be able to determine the time response 
of the chosen thermocouple a 50 mm x 50mm 
square shock tube arrangement was used. The 
shock tube consisted of two sections, a closed 
high-pressure driver section and an open ambient 
pressure driven section, with a thin plastic 
diaphragm separating the two sections. By filling 
the driver section with high-pressure air until the 
plastic diaphragm bursts a shock, or blast, wave 
is produced that may travel at speeds greater than 
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Table 2 Polynomial coefficients for ITS-90 
standard. 


Type K 
polynomial Value 
coefficients 
i d; 
0 0.226584602 
1 24152.10900 
2 67233.4248 
3 2210340.682 
4 —860963914.9 
5 4.83506 x 10'° 
6 —1.18452 x 10!" 
7 1.38690 x 10% 
8 —6.33708 x 101° 


the speed of sound (approximately 330ms'! 


in air at standard temperature and pressure) 
through the driven section of the tube. At the exit 
of the driven section of the shock tube a small 
heat source (flame) is positioned with the test 
thermocouple situated close to its core such that it 
measures the local flame temperature. When the 
shock wave exits the driven section the flame is 
extinguished and the thermocouple experiences a 
step change in temperature from that of the flame 
(approximately 466 °C) to that of the ambient 
temperature (approximately 21°C), Figure 4. 
The response time of the thermocouple is the 
time taken for the temperature to drop from the 
maximum temperature (7,,,,) to 63.2% of the final 


600 
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feos Tmax 
400 


Temperature (°c ) 
wo 
=) 
—) 


0 1000 2000 3000 4000 5000 6000 7000 
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temperature. Several 0.0785 mm (40-gauge) K- 
type thermocouples were evaluated and their time 
constants determined to be approximately 10 ms. 


Sparkler temperature profile analysis 


Initially, tests were performed on commercially 
available sparklers that had a nominal diameter 
of 2.8 mm and prepared to a length of 40 mm, 
and located horizontally in the apparatus shown 
in Figure 1. Since the surface temperature was 
known to be as high as 1575 °C, as measured 
previously with an R-type thermocouple, the 
locations of the K-type thermocouple were 
examined and it was concluded that they should 
be positioned no closer than 1.2mm from the 
sparkler surface, and a typical temporal profile 
is shown in Figure 5. Here, it may be observed 
that as the flame approaches the thermocouple 
bead there is a steady rise in temperature until 
the flame is adjacent to the thermocouple where 
the maximum temperature is recorded. Once the 
flame passes the thermocouple the temperature 
falls to ambient conditions. However, due to the 
instability and pulsating nature of the flame an 
oscillation in the temperature profile can also be 
noted in this Figure. 


Wasmann*® previously observed this fluctuating 
behavior of a pyrotechnic system and provided 
the explanation that two major factors were 
responsible for this action: 


(1) The competition between the various chemical 
reactions within the pyrotechnic composition. 


(2) Physical factors like heat loss, heat 
accumulation, and the intermittent vaporization 


” Wasa 


Local Temperature ( °C) 


Time (seconds) 


Figure 4 Thermal time response for a 40-gauge 
K-type thermocouple. 
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Figure 5 Zemperature profile taken 1.2 mm away 
from a sparkler surface. 
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processes in the test composition. 


In the first case, the pulsating flame behavior 
occurs as a result of competition in the chemical 
reaction and this may be explained as oscillations 
between light and dark cycles that occur when 
the composition consists of two or more different 
fuels. The quick reacting fuel reacts with the 
oxidizer releasing little energy in the form of 
light, hence representing a relatively dark cycle, 
whereas the light cycle is caused by large energy 
releases from the oxidation of the slow reacting 
fuel thereby producing large quantities of visible 
light. The rate of oxidation of the slow reacting 
component increases rapidly after the dark cycle 
due to the increase in the reaction surface area, 
and also due to the presence of oxidative gases 
from the dark cycle, trapped in the micro-porous 
structure of the surface. This dual cycle process 
seems to cause the pulsations in reactions with 
multiple fuel components, and with variable rates 
of reaction. In the case of sparklers, clear dark 
and light zones were not observed; however, the 
intensity of the emitted light could be seen to 
vary, and this pulsating behavior is observed as a 
change in emitted light intensity. In the chemical 
composition of the sparkler, the effect of the minor 
components has not been considered although 
these minor components may also affect the 
reaction rate by reacting with the oxidizer, thereby 
resulting in the dark phase of the pulsation cycle. 


In the second case, the undulating nature of the 
temperature profile may be explained by the rapid 
heat loss from the surface following the release of 
an iron particle. This rapid thermal change would 
be sufficient to instantaneously lower the local 
temperature. Once the particle was released the 
surface temperature would again return to a quasi- 
steady state where the local temperature would 
attempt to return to its original value, albeit lower 
in value because of the passing of the flame front, 
before the release of another iron particle. The rate 
of release of the iron particles would provide the 
fluctuations in the temperature profile. 


Although other physical properties like viscosity 
and volatility can also affect the burning process 
to create a pulsation effect as explained by 
Gol’binder and Goryachev,’ they were not 
considered in this study since there are no volatile 
components involved in the formation of the 
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standard sparkler. 


In order to determine the temperature distribution 
around the surface of the sparkler, measurements 
were taken of the maximum attainable temperature 
at discrete locations in the +x, and +z and —z 
directions, where +x refers to the horizontal 
distance radially away from the sparkler surface, 
and the +z and —z directions refer to the upper 
and lower vertical distances away from the 
sparkler surface respectively. Figure 6 shows 
diagrammatically the coordinate system adopted 
for these sets of measurements. These maximum 
temperatures were obtained in the three directions 
from more than 120 individual tests for a range of 
positions, from 1.15 mm to 2.5 mm in increments 
of 0.05 mm in the +x direction and 0.10 mm in the 
+z and —z directions, Figures 7 and 8. 


Although there is a large amount of scatter in the 
results, mainly due to the small differences in 
the composition of the charges, and the location 
of the thermocouple in relation to the sparkler 
surface, there is a distinct trend in the temperature 
distribution. From Figures 7 and 8, it may be 
observed that a drop in temperature occurs between 
1.2 mm to 2.5 mm in both the x and z directions. 
Comparing the curve-fitted temperature profiles 
in the x and z directions, it may also be observed 
that the drop in the temperature in the z (vertical) 


directions is estimated to be some 170°C mm‘, 


Mounting 
plate 


Figure 6 Coordinate system for the measurement 
of temperatures from the main charge (sparkler). 


Page 15 


whereas in the x (horizontal) direction the fall 
in temperature is approximately 3 times greater. 
Furthermore, the temperature profiles shown in 
Figure 8 suggest that the thermal gradient in the 
+z direction is less than that in the —z direction, 
and may be explained by the following: 


(1) As the pyrotechnic materials melt in the 
reaction zone they tend to flow producing a 
downward shift in the sparkler’s position due 
to its weight before complete combustion can 
occur. This effect reduces the distance between 
the thermocouple and the reacting pyrotechnic 
surface becoming more pronounced when the 
linear burning speed of the reaction zone, as 
it moves along the test sample, is low. In the 
current condition a burn speed of 2 mms ! 
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Figure 7 Temperature distribution in the horizontal + x 
direction. 
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Figure 8 Maximum temperatures above and 
below the sparkler surface. 
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was recorded. 


(2) Effect of natural convection on the thermal 
field; similar to that observed around a heated 
horizontal cylinder by a Schlieren technique. 


One of the difficulties of using standard sparklers 
was that of igniting the specimen accurately at a 
preset time. To overcome this obstacle, an electric 
match was formulated from two pyrotechnic 
components and attached to a main charge 
fabricated from sparkler material. 


Test charge composition 


The pyrotechnic material used in the main charge 
was obtained from commercial sparklers without 
the addition of the iron particles. The reason 
for removing these particles was to reduce the 
fluctuations in the temperature measurements 
previously observed and shown in Figures 3 and 5. 
This was achieved by crushing the sparklers into a 
fine powder and removing the iron particles with an 
electric magnet. The powder was then mixed with 
a chemical binder (dextrin) and moulded into a test 
charge, referred to as a commercial pyrotechnic 
material (CPM), The composition of this charge is 
given in Table 3: barium nitrate makes up 74% of 
the charge, and aluminum 24%, with 2% taken up 
with a chemical binder (dextrin). 


Since it was intended to ignite the test charge 
in a more systematic manner, an electric match 
was fabricated with the main test charge. The 
pyrotechnic composition of the electric match 
is shown in Table 4: the oxidizer was potassium 
chlorate and the fuel was powdered lead 
thiocyanate. To hold the pyrotechnic compositions 
together a binder (dextrin) was again utilized. 


Test charge fabrication 


Once the main charge comprising the barium 
nitrate, aluminum and its binder was thoroughly 
mixed, it was pressed into a small cylindrical 
mould that had an internal diameter of 4mm 
diameter, and was 15 mm long, for a length of 
12 mm, with the remaining 3 mm allocated to the 
material for the ‘electric match’. A steel support 
bar, 60 mm in length and 2 mm in diameter was 
then pushed through the mixture along the axis of 
the mold. An electric heating element consisting 
of a 50.8 mm long, 0.16 mm diameter (34-gauge) 
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Table 3 Chemical formulation of commercial pyrotechnic material. 


Charge Chemical Formula Quantity (by wt) Mesh size 
: Barium nitrate Ba(NO;), 74% <200 
Main charge 
Aluminum Al 24% <200 
(CPM) . . 
Binder (dextrin) (C6H) 90s), 2% <200 
Table 4 Chemical composition of electric match. 
Charge Chemical Formula Quantity (by wt) Mesh size 
Potassium chlorate KCIO; 55% <200 
Electric match — Lead thiocyanate Pb(SCN), 44% <200 
Binder (dextrin) (C¢6H,90s), 1% <200 
nichrome wire with a resistance of 2.73 ohm was Testing facility 


wrapped around the support bar but isolated from 
it by a small piece of electrical insulating tape. The 
electric match mixture was then pushed into the 
remaining 3 mm part of the mold encompassing 
the electrical element such that the flat surface 
of the match composition is in contact with the 
composition of the main charge, around which the 
temperature distribution was to be measured. The 
complete test charge was then allowed to dry for 
24 h and then removed from the mold. Ignition of 
the test charge was accomplished by igniting the 
electric match using the heating element connected 
toa5 V DC power supply. Figure 9 is an image of 
the complete test charge showing the main charge, 
electric match, support bar and the electrical leads 
of the heating element. 


Figure 9 /mage of fabricated CPM charge from 
sparkler material. 
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In order to improve the ability to locate a 
thermocouple bead away from the surface of the 
test charge with more accuracy a 3-dimensional 
traversing mechanism was constructed, 
Figure 10. 


The base of the main traverse section supported 
the test charge holder, test charge, and three 
micrometer controlledslides to which three separate 
thermocouples could be mounted. The slides 
were configured to provide accurate independent 
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Figure 10 Facility for accurately positioning 
thermocouples close to the main charge. 
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movement in three orthogonal directions (x, y, 
z), with a spatial resolution in each direction of 
1 micron. To allow the test charge to be rigidly 
fixed in space a separate aluminum mounting 
plate was utilized, and two extension arms added 
ensuring that the test charge was mounted off 
the plate by approximately 40 mm, Figure 10. In 
addition, a second flat rectangular traversing block 
was attached to the vertical segment of the traverse 
and acted as the base for multiple thermocouple 
attachment, as shown in Figure 10, thereby 
allowing multiple measurements to be made for 
any one test. 


CPM temperature profile measurements 


With a test charge mounted in the extension arms 
of the mounting plate an electrical connection to a 
5 V DC power supply can be made to initiate the 
ignition process. Temperature measurements were 
taken using three K-type thermocouples arranged 
so that measurements could be obtained in the x- 
direction, and the +z directions simultaneously 
at precise horizontal and vertical locations, 
Figure 6. 


In this phase of the study the flame of the 
pyrotechnic was larger than the 2.8 mm diameter 
sparkler, which resulted in the burnout of the K- 
type thermocouples at approximately 2 mm from 
the surface of the test piece. Therefore in order to 
prevent continual thermocouple failure a larger 
standoff distance was utilized when compared 
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Figure 11 Temperature distributions around a 4 
mm diameter CPM charge. 
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to the 2.8mm diameter test charge. Figure 11 
shows three temperature profiles around a 4mm 
diameter test charge and it can again be observed 
that nearer the surface of the test charge, the 
magnitude of the temperature profile beneath the 
burning pyrotechnic (in the —z direction) is larger 
than that above the charge (in the +z direction), a 
trend similar to that found for the 2.8 mm diameter 
test species. However, the most striking feature 
of this data set is that the drop in temperature 
away from the burning test charge falls in a tight 
band with a nominal 200 °C spread for all three 
directions, providing a temperature decay of 
approximately 240°C mm‘! over the distances 
measured. This is in contrast to the results found 
for the 2.8 mm diameter sparkler. However, there 
are small differences between the data sets and it 
may be inferred that the maximum decreases in 
temperature for all the three orientations occur 
in the x directions. The minimum drop may be 
observed in the +z direction and the reason may 
be attributed to the flow of gases, due to natural 
convection, in this direction. The temperature drop 
in the —z direction is between the temperature drop 
values in the other two directions, but below that 
for the +z direction. The reasons for this relatively 
lower temperature drop in the —z direction can be 
attributed to interference with natural convection 
from jet-like ejections, not related to the emission 
of iron particles. 


The primary reason for the scatter in the data 
was due to the slight uncertainty in thermocouple 
position and the continuous pulsations in the flame 
when the premixed composition burned. Other 
factors that influence data scatter are related to the 
following initial assumptions: 


e Each test piece has the same density/ 
composition. 


¢ Burning is uniform — following a ring-like 
pattern along the length of the test piece. 


e Each test piece is unaffected by external 
conditions. 


Although care was taken to maintain constant 
weight of each test charge, their weights varied 
within +2% error. Similarly while compacting the 
composition into the mold shell, attempts were 
made to maintain the dimensions to within +2%, 
thereby keeping the test charge final density to 
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Table 5 Heat of formation for different chemical products. 


Chemical formula Chemical name 


Ba(NO3), Barium nitrate 
Al Aluminum 

BaO Barium oxide 
ALO; Aluminum oxide 
N, Nitrogen gas 


within a couple of percent. 


Furthermore, it was assumed that the ignition of 
the main charge occurred instantaneously from 
the initiator (electric match) with the resulting 
flame propagating uniformly through the entire 
cross-section of the main charge. In the absence of 
pulsations a ring-like burning pattern was assumed. 
However, if the flame does not cover the entire 
cross-section, an irregular burning pattern results 
and this can be exacerbated by the difference in 
the surface burning speeds. It was also assumed 
that changes in the ambient conditions, such as 
variations in temperature and forced convection, 
were small enough to be neglected. 


It should also be noted that no temperature 
corrections due to conduction, convection and 
radiation have been applied to these measurements 
since it was considered that any error analysis 
would have been dependent on the thermal 
conductivity of the wires (A), the convective 
heat transfer coefficient (4) and the emissivity 
(e) of the thermocouple sensor element. This 
dependence would have relied upon a steady state 
energy balance, a condition that was inappropriate 
for this dynamic situation. Furthermore, due to 
the gaseous nature of the pyrotechnic event, the 
convective heat transfer coefficient could not have 
been accurately determined; neither would the 
emissivity of the thermocouple bead be a single 
value due to the heating and cooling nature of the 
event discolouring the bead surface. However, as 
an approximation, a conduction error has been 
estimated from a simple energy balance between 
the net heat conducted along the wires and the 
heat convected to the wires. This analysis forms a 
simple second order differential equation that can 
be solved for the sensor temperature. For example, 
for a thermal conductivity of 24Wwm'' deg! 
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Heat of formation AH; /kJ mol! 
=997 
0 
—548.1 
—1675.7 
0 


for the wires, and 45.4Wm'deg'! for the 
convective heat transfer coefficient, an error 
in the measurement of a bead temperature of 
1200 °C would have been approximately 130 °C, 
an error of some 10.8%. Likewise, an error in 
the temperature of the bead due to radiation has 
been estimated from a simple steady state energy 
balance between the convection heat transfer and 
the radiation transfer. This was achieved by using 
the above value for the convection heat transfer 
coefficient, and an emissivity value of 0.07 for a 
K-type (chromel/alumel) thermocouple, where it 
was found that the error in the temperature was 
approximately 1860 C, or 15.5%. 


Summary 


This work has shown that the near surface 
temperatures of a pyrotechnic charge may be 
measured accurately with miniature high-speed 
K-type thermocouples. Temperatures approaching 
1600 °C were recorded at the surface of a 2.8 mm 
sparkler using an R-type thermocouple, while 
temperatures over 1300 °C were being measured 
1.2 mm away from the sparkler surface using K- 
type thermocouples. 


In the case of a standard sparkler the temperature 
distribution may be measured as close as 1.2 mm to 
the surface, whereas for the 4 mm diameter charge, 
made from the same material, the thermocouples 
could not be placed within 2.0 mm of the surface 
without a total failure. 


When comparing Figures 7 and 8 for the 2.8 mm 
diameter sparkler to Figure 11 for the 4 mm charge, 
the temperature decay for the sparkler appears to 
have directional dependence and falls off more 
slowly than that for the 4mm charge. It was also 
found that for the 2.8 mm sparkler the temperature 
decay in the vertical direction was approximately 
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30% lower compared to the temperature decay 
in the horizontal direction. This suggests that the 
heat liberation rate for the 2.8mm charge was 
insufficient to overcome the energy transport 
by natural convection and thereby resulted in 
lower levels of thermal energy propagated in the 
horizontal direction. 


The temperature around a 4mm diameter test 
charge fabricated from commercially available 
pyrotechnic material has shown that, although 
there is a difference in the thermal energy being 
released in different directions, this difference 
is small compared with that obtained from the 
2.8 mm sparkler, and that the temperature falls off 
in a tight band of approximately 200 °C. 
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Abstract: An accidental explosion in 1848 in the USA became a trigger for the development of neuroscience. 
An accidental explosion of black powder took place in a borehole for blasting and the expelled iron bar 
penetrated the head of a young man. He was injured but survived for 12 years. The authors were asked to 
calculate the speed, impact pressure or energy of the explosion by the producer of a TV program. At the time 
we were carrying out similar experiments using a mortar and firework stars, and so a model experiment 


was performed. Here we report of the results. 


Keywords. black powder, ballistics, iron bar, shot energy efficiency 


Introduction 


An accidental explosion involving — black 
powder became a trigger for the development 
of neuroscience. On September 13, 1848, at a 
blasting site in Vermont, USA, Phineas Gage, a 
young man of 25 years old, was involved in an 
accident.' He accidentally dropped an iron tamping 
bar 6.1 kg in mass, 3.2 cm in diameter and 110 cm 
in length into a borehole 3.8 cm in diameter and 
91 cm deep filled with about 0.4 kg black powder. 
An explosion took place, and the bar was shot 
from the borehole and penetrated his head from 
his left cheek to the top of his head. Amazingly he 
escaped death and survived for twelve years after 
the accident. 


Twenty years after the accident, Gage’s physician, 
John Harlow, perceptively correlated Gage’s 
cognitive and behavioral changes with a presumed 
area of focal damage in the frontal region. His 
observation made considerable scientific impact 
and gave rise to controversy. Gage’s skull was 
recovered and is in the Warren Anatomical Medical 
Museum at Harvard University. 


We were asked by the producer of a television 
program to calculate the velocity of the iron bar, 
the pressure of the bar on his brain, and the power 
of the explosion. In response to his request we 
carried out a small scale experiment using a mortar 
for firework star shots and a pressure measuring ap- 
paratus? and the results were analyzed as before.* 
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Experimental 
Materials 


The grain black powder and electric match were 
made by Nippon Kayaku Co. Ltd. This black 
powder may be similar to the black powder used in 
blasting 160 years ago in the USA. Howard? shows 
that the black powder manufacturing process has 
not changed in principle from ca. 1780 until now. 
Van Gelder and Schlatter® wrote in 1927 that 
blasting powder was made in very much the same 
way as gunpowder. 


Apparatus 


The iron bar is 1.5 kg in mass, 22 mm in diameter 
and 506 mm long. The mortar is shown in Figure 1. 
The inner diameter, wall thickness and depth of 
the mortar were 25mm, 2.5mm and 460 mm, 
respectively. The gap ratio in this case was 29%. 
The gap ratio GR (%) is defined as follows: 


Here, S and s are the cross sectional areas of the 


S-s 


GR= x 100 


mortar and the bar, respectively. 


Two pressure sensors were fitted, to the bottom 
and to a position 300 mm from the bottom of the 
mortar. The mortar was fixed on an H-shaped steel 
holder as shown in Figure 2. 
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~~ 12V battery 
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Figure 1 Mortar and bar. The pressure in the mortar during the shot was measured and recorded using 
two pressure sensors (Kistler 60410A), charge amplifiers (Kistler 5011) and a digital oscilloscope (Sony 


Tektronix 5011). 


Procedure 


The mortar was set on the ground vertically. The 
black powder and electric match were put in the 
bottom of the mortar. Then, a piece of tissue paper 
and the iron bar were inserted into the mortar 
slowly. The mortar fixed to the holder was set on 
the ground horizontally. The mortar was covered 
by U-shaped concrete blocks for safety. The space 
in front of the mortar was protected by sand bags 
and concrete blocks as shown in Figure 2. 


The electric match was ignited by an electric 
current. The black powder burned, pressure 
developed and the bar moved forward. The 
pressure profile was recorded on an oscilloscope 
and the muzzle velocity of the bar was estimated 
from the pressure profile. 


Mortar afid iron 
bar 


Sand bag 


Figure 2 Setup of the shot experiment of an iron bar. 
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Results and Discussion 
Pressure profile 


An example of the observed pressure profiles in 
the mortar while the bar was shot is shown in 
Figure 3. The pressure profile shown by the fine 
line was recorded by the bottom sensor, and the 
profile shown by the thick line by the middle 
sensor. Both profiles in the figure were recorded 
simultaneously in the shot. 


As indicated later, the bar stayed in the mortar while 
the pressure developed, that is, the combustion 
of the black powder had finished before the bar 
left the muzzle of the mortar. The pressure of the 
bottom sensor is higher than that of the middle 
sensor. The pressure of the bottom sensor is that 
in the space of the mortar behind the bar, and the 
pressure of the middle sensor is that in the gap 


Pressure measuring 
apparatus and ignition 
divice 


U- shape 


H- shape steel concrete block 


holder 
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14 7 —at bottom 
—— 300mm from bottom 


Pressue(MPa) 
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Figure 3 Pressure profiles in the mortar with 
7.5 g black powder. 


between the mortar wall and the bar. 


Motion of a bar in the mortar 


The motion equations of a bar in the mortar are 
expressed as follows: 


du 

Se As 
ar p(t) (1) 
dz — 
eo (2) 


Here, m, u, A and z are the mass, velocity, cross 
sectional area and traveling distance of the bar, 
respectively. 


as (3) 


Here, D is the diameter of the bar, and 


fu 2D” nl) 
dt 4m . (4) 


Here, p(t) is the observed value and is substituted 


into equation (4) 


Equations (2) and (4) are simultaneously solved 
by numerical calculation, and acceleration du/d¢, 
velocity u and traveling distance z are obtained. 


Equations (2) and (4) were solved by the Runge— 
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Kutta method. The time integration process for 
ordinary differential equations (2) and (4) was 
performed using a fourth order accuracy Runge— 
Kutta method. 


The digital pressure data were recorded on 
an oscilloscope and the data were reduced 
using Excel. These reduced data were used for 
calculating acceleration, velocity and the distance 
traveled by the bar. 


The calculated profiles of the acceleration, velocity 
and distance of the bar are shown in Figure 4. 


In all cases, the velocity of the bar increased with 
time at first, then became constant, suggesting that 
the combustion of the black powder was complete 
and that the overpressure in the mortar reached 
zero. This situation can also be seen in time— 
distance curves in Figure 4. The distance of the bar 
from the bottom of the mortar increases with time, 
but at a point, the rate of increase changes, and 
then it becomes lower and constant. The change 
point occurs when the inner overpressure in the 
mortar reaches zero. With 10 g black powder, the 
inner overpressure reached zero just before the bar 
left the muzzle. 


The observed and calculated results of the bar shot 
experiment are listed in Table 1. 


Effect of mass of black powder on peak 
pressure and muzzle velocity 


The effect of the mass of black powder on the peak 
pressures in the mortar is shown in Figure 5. The 
peak pressures at both the bottom and the middle 
of the mortar increased exponentially with the 
mass of black powder. 


Plots of the calculated muzzle velocity and kinetic 
energy of the bar against the mass of black powder 
are shown in Figure 6. The calculated muzzle 
velocity increased linearly with the mass of black 
powder (BP), but the calculated kinetic energy 
showed an exponential increase with the mass of 
BP. 


Energy efficiency of the shot of the iron bar, 
shell and star from mortar 


It is necessary to know the energy efficiency of the 
shot of the iron bar from the borehole in Gage’s 
case, in order to estimate the initial velocity and 
kinetic energy of the bar. The shot energy efficiency 
of the iron bar, shell and star which were observed 
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Figure 4 Calculated profiles of the velocity and distance of the bar in the mortar. 


in experiments are listed in Tables | and 2. The 
shot energy efficiency in Tables 1 and 2 is the 
kinetic energy of the shot of a projectile divided 
by the explosion energy of black powder. The 
explosion energy of black powder was measured 
and published as 2800 J g ' by Rose.’ 


We could not perform the experiment on the same 
scale as Gage’s case, because no safe facility was 
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available for the experiment. The information from 
this experiment was limited so it was necessary to 
estimate the data in Gage’s case. So we used for a 
supplement the star and shell data which we had 
previously obtained.** 


A plot of the shot energy efficiency of the iron 
bar against the mass of BP in this work is shown 
in Figure 7. In this case, with 23% GR, the shot 
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Table 1 Observed and calculated results (GR = 23%). 


MassofBP Praxis Zuma mazda AT CMV. KE EE 
Run (g) (MPa) (ms) — (MPa) _(ms) (ms) (ms') (J) (%) 

1 15 = = 

5 10 23.3 2.88 11.86 4.06 1853 248 460 1.64 

3 7.5 13.55 5.54 6.94 6.86 20.84 17.4 227 1.08 

4 5 5.42 7.86 77 8.74 24.1 10.9 89 0.64 

5 2.5 2.4 9.22 1.16 10.73 25.66 5.4 ye) 0.31 


BP: black powder; P,,,.;: Maximum pressure by the bottom pressure transducer; P,,,..9: Maximum pressure by 
the middle pressure transducer; T,,,,: time to maximum pressure: AT: time during positive overpressure; CMV: 
calculated muzzle velocity; KE: kinetic energy of the bar; EE: energy efficiency of the shot 
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Figure 8 Plot of the shot energy efficiency of a 
no. 3 firework shell vs. BP mass with about 13% 
BP mass(g) GR. 


Figure 6 Plots of the calculated muzzle velocity 
and kinetic energy of a bar. 
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Table 2 Shot energy efficiency of iron bar, shell and star. 


(1) No. 3 Shell. Dimensions of mortar: @ 90 mm x 750 mm (Ref. 4,9) 


Muzzle 

GR velocity KE EE 

(%) (ms ') (J) (%) 

13 135 2208 2.63 
12 122 1797 2.29 
13 132 2061 3.07 
12 82 807 1.31 
14 78 730 1.30 
13 89 958 1.90 
13 83 834 1.86 
14 88 926 2.36 
14 70 581 1.73 
15 73 679 2.43 
15 57 419 1.87 
15 30 113 0.67 
15 23 66 0.39 
13 33 131 0.78 
15 34 139 0.83 


(2) Star. Dimensions of mortar: @ 25 mm = 460 mm (Ref. 8) 


Shell 
BP mass Shellmass diameter 
Run (g) (kg) (m) 
1 30 0.242 0.084 
2 28 0.242 0.084 
3 24 0.237 0.084 
4 22 0.240 0.084 
5 20 0.240 0.084 
6 18 0.242 0.084 
7 16 0.242 0.084 
8 14 0.239 0.084 
9 12 0.237 0.084 
10 10 0.255 0.083 
11 8 0.258 0.083 
12 6 0.250 0.083 
13 6 0.250 0.083 
14 6 0.241 0.084 
15 6 0.241 0.083 
BP mass Star mass Star diameter 
Run (g) (kg) (m) 
1 1 0.0095 0.0224 
2 0.4 0.0095 0.0227 
3 0.7 0.0091 0.0223 
4 2 0.0092 0.0222 
5 1.5 0.0083 0.0218 
6 1.5 0.0089 0.0224 
7 2 0.0086 0.0225 
8 1.5 0.0088 0.0226 
9 1.25 0.0083 0.0224 
10 1.75 0.0085 0.0220 


energy efficiency increased linearly with mass of 
BP. 


A plot of the shot energy efficiency of a no. 3 
firework shell against BP mass is shown in Figure 8. 
The shot efficiency increased with BP mass, 
though the observed data are very scattered. 


Figure 9 shows a plot of the shot energy efficiency 
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Muzzle 

GR velocity KE EE 

(%) (m/s) (J) (%) 
20 87 36 1.29 
18 27 3 0.31 
20 57 15 0.75 
21 117 63 1.13 
24 103 44 1.05 
20 106 50 1.19 
19 117 59 1.05 
18 86 33 0.78 
20 77 25 0.70 
22 91 35 0.72 


of firework stars. With about 20% GR, the shot 
energy efficiency increased with BP mass in the 
range of 0-2 g BP. But with about 30% GR, the 
efficiency decreased somewhat with increasing 
BP mass when the BP mass exceeded 2 g. 


Figure 10 shows a plot of the shot energy efficiency 
of stars vs. GR with 2 g BP mass. The efficiency 
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Table 2 continued Shot energy efficiency of iron bar, shell and star. 


(3) Star. Dimensions of mortar: © 20 mm x 360 mm (Ref. 8) 


BP mass Star mass Star diameter 


Run (g) (kg) 
1 1 0.0038 
2 a 0.0038 
3 3 0.0037 
4 4 0.0041 
5 5 0.0033 
6 6 0.0038 
7 7 0.0037 
8 2 0.0039 
9 2 0.0040 
10 2 0.0035 
11 2 0.0041 
12 2 0.0038 
13 2 0.0019 
14 2 0.0018 
15 2 0.0018 
16 2 0.0020 
17 2 0.0018 
18 2 0.0009 
19 2 0.0008 
20 2 0.0009 
21 2 0.0011 
22 2 0.0010 


(m) 

0.0167 
0.0172 
0.0169 
0.0174 
0.0166 
0.0167 
0.0171 
0.0171 
0.0175 
0.0162 
0.0175 
0.0167 
0.0134 
0.0133 
0.0132 
0.0138 
0.0132 
0.0103 
0.0105 
0.0104 
0.0111 
0.0110 


© GR 20% 
OGR 30% 
AGR 55% 
X GR 70% 


Shot energy efficiency (% 


BP mass(g) 


Figure 9 Plot of the shot energy efficiency of 
firework stars vs. BP mass with various GR. 
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Figure 10 Plot of the shot energy efficiency of 
firework stars vs. GR with 2 g BP mass. 
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Table 3 The maximum shot energy efficiency of iron bars, no. 3 firework shell and stars. 


Shell Star 
Projectile mass (kg) 0.24 0.0038 
BP mass (kg) 0.024 0.0020 
Projectile/BP 10 1.9 
GR (%) 13 30 
Shot energy efficiency (%) 3.1 0.54* 


Bar Bar 
Star (this work) (Gage’s case) 
0.0038 1.5 6.1 
0.0060 0.010 0.40 
0.63 150 15 
31 23 29 
0.41** 1.6 ? 


*Maximum shot energy efficiency with 2.0 g BP mass. **Not maximum shot energy efficiency with 6.0 g BP mass. 


decreases with increasing gap ratio. 


The information about the maximum shot energy 
efficiencies of iron bars, no. 3 shells and various 
stars are listed in Table 3. In this Table, the 
maximum shot energy efficiency was obtained 
with a small gap ratio (13-20%) and the largest 
BP mass in the experimental range. However, 
with a larger gap ratio (about 30%), the efficiency 
decreased with increasing BP mass (<2 g). 


Estimation of the velocity and kinetic energy 
of the iron bar in Gage’s case 


At the moment it is difficult to estimate the shot 
energy efficiency in Gage’s case accurately. Most 
of data showed that the efficiency increased with 
BP mass in the experimental range. But with 
about 30% gap ratio the efficiency decreased 
with BP mass when the mass exceeded a critical 
value. From above considerations, the shot energy 
efficiency of the iron bar in Gage’s case may be 
roughly estimated at about 1%. 


The kinetic energy and muzzle velocity of the 
iron bar in Gage’s case are estimated to have been 
about 11 kJ and 60 ms |, respectively. A more 
reliable estimation may be obtained by conducting 
a similar experiment using the same experimental 
conditions as in Gage’s case. 


In our experiment with 15 g black powder, the 
expelled iron bar penetrated a 20cm long sand 
bag and then made a hole in a concrete board 
about 5 cm thick. 


According to the Department of Defense, USA,'°a 
hazardous fragment is one having an impact energy 
of 79 J or greater. The energy of the expelled bar 
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in Gage’s case was much greater than the energy 
of a hazardous fragment. 
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Ballistics of a No. 3 Spherical Shell with Illuminant 
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Abstract: Shot experiments were carried out using a No. 3 spherical shell with an illuminant. The three 
dimensional trajectory of the flying shell was obtained by tracing the shell with two high-speed video 
cameras in different orientations and by analyzing the recorded video picture. The 3DSTARI computer 
code was developed for calculating the three dimensional trajectory of a flying shell with high accuracy. 
The optimal drag coefficient, wind speed and wind direction were estimated using the 3DSTARI code by 
fitting the calculated trajectory to the experimental one for a No. 3 shell. 


Keywords: ballistics, no. 3 shell, trajectory 


Introduction 


A firework shell is shot from a mortar, rises into 
the air, bursts in the sky and releases burning stars. 
The stars fly into the sky and form a composition 
like a blooming flower or bunch of flowers. The 
shell is set in the mortar on top of a lifting charge. 
The lifting charge is ignited and the shell is shot 
into the air by the pressure developed by the 
combustion gas of the lifting charge. 


The ballistics of the shell expelled from the muzzle 
of the mortar are affected by various factors such 
as muzzle velocity, air drag, shot angle, mass of 
lifting charge, wind direction and speed, and so 
on. In terms of the design and safety of the shell 
shot, it is important to know the basic ballistics of 
the shell. 


Shimizu!” has carried out shot experiments using 
spherical shells and analyzed the results. Kosanke 
and Kosanke* have performed _ theoretical 
modeling of the ballistics of shells using Shimizu’s 
experimental data. Recently, lida et al.’ have 
carried out a shot experiment using several sizes 
of spherical shells. Eckhardt and Andre,°> and 
Speer® have calculated the trajectories of spherical 
firework shells in order to investigate an accident 
at a public fireworks display in 1997. Mercer’ has 
modeled the aerodynamics of propelled aerial 
shells. Schneider and Schneider* performed 
ballistic trajectory calculations to investigate 
the relationship between the launch elevation of 
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dud aerial firework devices and ground impact 
distances. 


The objectives of this paper are as follows: 


(1) A No. 3 spherical shell is shot from a mortar 
and the trajectory of the shells is observed 
from different directions by two high-speed 
video cameras, and the results are three- 
dimensionally analyzed. 


(2) A three-dimensional theoretical model of the 
ballistics of the shell is developed, a theoretical 
calculation is carried out, and the effect of 
wind direction and speed, air drag, shot angle 
and so on are examined. 


Experimental 
Materials 


No. 3 spherical shells with an illuminant for tracing 
the trajectory were made by the Sunaga Fireworks 
Company. The lifting charge and electric match 
were made by the Nippon Kayaku Company. The 
mass and diameter of the shells and the mass of 
lifting charges in this experiment are listed in 
Table | along with the observed muzzle velocity. 


Apparatus 


The mortar for the No. 3 shells was made of steel 
and the dimensions of the mortar are shown in 
Figure 1. 


The trajectory of an expelled shell was traced 
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Table 1 Mass and diameter of the shell, mass of lifting charge, and muzzle velocity. 


Run Mass Diameter 

(kg) (m) 

1 0.255 0.083 

2 0.258 0.083 

3 0.250 0.083 

4 0.250 0.083 

5 0.241 0.084 

6 0.241 


Lifting charge Muzzle velocity 
(g) (ms ') 
10 bie 
57 
30 
23 
33 


DNDN DW 


34 


0.083 


by two high-speed video cameras (Phantom VR- 
V4.2) with a frame rate of 150 frames per second. 
The locations of the two cameras are shown in 
Figure 2. 


Procedure 


The mortar was set on the ground vertically. The 
lifting charge and electric match wrapped in thin 
paper were put in the bottom of the mortar. Then 
a No. 3 shell was placed on the lifting charge. 
The electric match was ignited by turning on an 
electric current. The lifting charge burned, pressure 


09m, 


0.78r 


Mortar 


Lifting charge with 
electric match 


Figure 1 Dimensions of the mortar for No. 3 
shells. 
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developed, the shell moved upwards and was 
expelled into the air. The trajectory of the shell in 
the air was recorded by the high-speed cameras. 


Three dimensional analysis of experimental 
data 


The position of the shell flying in the air is 
analyzed using two cameras located as in 
Figure 2. Two three-dimensional rectangular 
coordinates are provided for the locations of the 
cameras. The muzzle of the mortar is the origin 
of the coordinates. The perpendicular direction is 
the z axis of the coordinates. The two horizontal 
coordinates are X—Y and x—y, and camera | is on 
the X axis of the X—Y coordinates and camera 2 
is on the y axis of the xy coordinates as shown 
in Figure 2. The Y coordinate of the shell in the 
X—Y coordinates is recorded by camera | and the 
x coordinate in the x—y coordinates is recorded by 
camera 2. The relationships of the two coordinates 


xX 


Camera 1 ———__ 


Camera 2 


206.95m 


Mortar 


139.23m 
Figure 2 Location of high-speed cameras. 
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are as follows. 


x= X cosa —Ysina (1) 


y=Xsina+Ycosa (2) 


The coordinate y is calculated from the angle a, 
and the recorded x and Y values, using equations 
(1) and (2). 

Finally, the spatial position (x, y, z) of the firework 
shell in one set of three-dimensional rectangular 
coordinates is obtained. These coordinates were 
converted to the real distances, and the real spatial 
position of the shell at a given time is calculated. 


Three-dimensional model for 
ballistics of a shell 


Equations of motion 


As a shell flies into the air, a complex aerodynamic 
drag force acts on it. The force will depend on the 
density of the air, the viscosity of the air, the shape 
and surface roughness of the shell, etc. To simplify 
the problem, the following assumptions are made. 


(1) The shell does not spin in flight. 
(2) The shell is spherical. 
(3) The air density and air viscosity are constant. 


(4) The speed and the direction of the wind are 
constant. 


The vectors of position, velocity and acceleration 
of the shell, and the force acting on the shell, 
are expressed by rectangular coordinates. For 
example, the vector of the position of a shell is 
expressed by rectangular coordinates as shown in 
Figure 3. 
F=xi +yj+zk 
x=rsinf, cosa, y=rsinB, sina, 
z=reos§, 
reax’+y 427 (3) 


The motion equation of a spherical shell is 
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Figure 3 Three dimensional coordinates of 
positional vector F , 


tana, =~ 
Xx, 


| 2 2 
tan B, Ne 


Z 


Here, i , r and k are the unit vectors of the 
coordinates x, y and z, respectively. 

a, is the angle between the shadow of the vector 
F in the x-y plane and the coordinate x, and £1 is 
the angle between the vector ~ and the coordinate 


Z 
expressed as follows: 


(4) 


Here, &, ¥, Fy and F’, are the vectors of 
acceleration under gravity, velocity of the shell, 
air drag and buoyancy, respectively. m is the mass 
of the shell. 


Furthermore, 
a dr 
di (5) 
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B (7) 


Here, p, p,, A and Cp are air density, density of 
the shell, cross-sectional area of the shell and drag 
coefficient of air, respectively. 7 and w are 
positional vector of the shell and relative velocity 
vector between air and the shell. 


The motion velocity of a shell: 


Vevitv,jtvyk (8) 


v, =|v|sin B, cosa, 
v= Iv sin B, sina, 
v, =|v\cosB, 

2 


2 2 2 
vo=v. tv, +y, 


x 


Vy 
tana, =— 
v 


x 


+v 


tanB, = - 
v 


z 
The wind velocity vector w is constant, 

w= w,t wy F wk (9) 
Here, 

w, =|nsin B, cosa, 

w, =|w{sin B, cosa, 


w, =|wsin B, cosa, 


Relative velocity between the shell and air is: 


(,-w i t+(v,-witQ,-w)k | 
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(v,-w,) +, -w,P+,-w,y OD 


Therefore, the motion equations of a flying shell 
in the air are: 


dv _ p 


)g—k|¥ —W|(¥-W) (12) 


In the three dimensional rectangular coordinate 
shown in Figure 3, the gravity acceleration can be 
expressed as follows: 


g =-gk 
And the equations of motion are: 

d ae (13) 

vs =—-k|v—w-(, —w,) 

dt _— 
dx (14) 
— — Vi 

dt 

dy (15) 
ra = kv - J (¥, = W,) 

(16) 

o, 

dt - 
dv Lo (17) 

2 =-(1- 2) gk | -i-(v, -w,) 

dt P, 
dz (18) 
—_— = Vv, 
dt 
Here, k = ie) ; when the shell is a sphere, 

2m 
3 
k= PE» , where D,, is the diameter of the 
40 yo? F 

sphere. 


Numerical calculation of motion equations 


It is difficult to integrate equations (13)-(18) 
directly. Therefore, numerical calculation of 
equations (13)-(18) was carried out using the 
fourth-order Runge—Kutta method in this work. 
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Results and Discussion 
Calculation accuracy of the 3DSTARI code 


The numerical calculation results by 3DSTARI 
were compared with the analytical solutions, and 
the accuracy of the former was validated. 


At first, considering the motion of the shell 
expelled perpendicularly into the sky without air 
drag, the perpendicular velocity V, and position 
z were calculated using 3DSTAR1. The muzzle 
velocity of the expelled shell was assumed to be 
119 ms | and the calculated results are shown in 
Figure 10. 


On the other hand, V, and z can be obtained 
analytically as follows: 


V.=Vo-gt (19) 


1 
za2ytVy-t- Sgt (20) 


As shown in Figure4, the results from the 
3DSTARI numerical calculation agreed com- 
pletely with the analytical solution of equations 
(19) and (20). 


Coodinate z {ml 


Effect of time interval on calculation accuracy 
of 3DSTARI code 

Kosanke and Kosanke* published ‘Computer 
Modeling of Aerial Shell Ballistics”, and validated 
the results using Shimizu’s experimental data. The 
effect of the time interval on the accuracy of the 
3DSTARI code was examined using the same 
data: 


Muzzle velocity of the shell: 118.95 ms! 
Diameter of the shell: 0.17125 m 

Mass of the shell: 2.1111 kg 

Angle of the mortar: 0° 

Wind speed: 0 ms ! 
Drag coefficient: 0.36 


When the drag coefficient is constant, the analytical 
solutions of the equations of motion of the shell 
expelled perpendicularly are as follows 


Time to 1 
apogee 


ty 


= ——. tan! Ley (21) 
gk g ° 


3 


8 


S 
Vebcity Vz(m/'s) 


= 


B 


Figure 4 Comparison of 3DSTARI calculation and analytical solution. Without air drag and with 


119 ms‘ muzzle velocity. 
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Apogee l k_, 

height Z,=—-Inj1+—V, (22) 
2k g 

Velocity 

a a & (ie? ) (23) 

impact k 


an nt ge ta( FE, (24) 
2 1 Jgk g 2 


Here, k = SPC and g, p, Cp, p, D, and Vy are 
ap 1). 


acceleration under gravity, air density, air drag 


coefficient, and density, diameter and muzzle 
velocity of the shell, respectively. The velocity of 
the shell is positive when the shell moves upward, 
and negative when it moves downward. 


The results are shown in Table 2 and Figure 5. 


It can be seen from Figure 5 that even if the time 
interval is 1 s, the computer code can give high 
enough calculation accuracy. 


Results of observed and calculated three 
dimensional trajectory 


The video pictures of the motion of the shell in 
the air were recorded by two cameras facing in 


Figure 5 Calculated results by 3DSTAR1. 


Table 2 Calculated results by 3DSTAR1 along with analytical solution. 


Time interval 


(s) (s) 
1.0 7 
0.1 7.0 
0.01 7.06 
0.001 7.063 
Analytical solution by equations 7.0637 


(30)-(33) 
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Time to apogee 


Apogee height Time to impact 
(m) (s) 
314.246 16 
314.3288 16.0 
314.3479 16.08 
314.3479 16.087 
314.3485 16.0857 
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Figure 6 Observed and calculated trajectories of 


an airborne No. 3 shell (run 1). 
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Figure 7 Observed and calculated trajectories of 
an airborne No. 3 shell (run 2). 


different directions. These pictures were analyzed 
and the trajectory of the shell was expressed using 
three dimensional rectangular coordinates. The 
results are shown in Figures 6-11. 


It is seen from Figures 6—11 that the shells were 
expelled nearly perpendicularly but moved with a 
tilt angle at high altitude. This may be the effect 
of the wind. 


The muzzle velocity of the shell was determined 
from the trajectory of the shell in the air. The 
results are listed in Table 1. The muzzle velocity 
tends to increase with increasing mass of lifting 
charge, but the scatter is large. 


Trial and error calculations were carried out using 
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Figure 8 Observed and calculated trajectories of 
an airborne No. 3 shell (run 3). 


® 8 


8 


Coordinates(m) 
a 


Figure 9 Observed and calculated trajectories of 
an airborne No. 3 shell (run 4). 


the 3DSTARI code for fitting the drag coefficient, 
wind speed and wind direction to the observed 
trajectory of the shell, and the optimal values were 
obtained. These values are listed in Table 3. The 
calculated trajectories with these values are also 
shown in Figures 6-11. 


Conclusions 


A shot experiment using a No. 3 spherical shell 
with an illuminant was carried out and three 
dimensional calculations were performed. The 
following conclusions were derived: 


(1) The three dimensional trajectory of the flying 
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Figure 10 Observed and calculated trajectories 


Figure 11 Observed and calculated trajectories 
of an airborne No. 3 shell (run 5). 


of an airborne No. 3 shell (run 6). 
shell was obtained by tracing the shell with 


two high-speed video cameras aimed in Acknowledgement 
different directions, and by analyzing the The authors wish to gratefully acknowledge 
recorded video picture three dimensionally. the experimental assistance of the Showarika 
(2) The 3DSTARI1 code was developed for Company, and re mh Aoyagi, an undergraduate 
calculating the three dimensional trajectory sua geue inte Pie aa baboraicry: 
of a flying shell with a high calculation 
accuracy. 
(3) The optimal drag coefficient, wind speed and 


wind direction were estimated using 3DSTAR1 
code by fitting the calculated trajectory to the 
experimental one for a No. 3 shell. 


Table 3 Muzzle velocity and drag coefficient of the shell, wind speed and wind direction. 


Muzzle velocit . Wind speed Angle between directions of 
is (ms ') , Dregcoemeient (m a a and coordinate x (°) 
1 73 0.53 5 225 
2 57 0.53 10 —110 
3 30 0.99 7 0 
4 23 0.71 4 10 
5 33 0.67 7 8 
6 34 0.63 8 —80 
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Abstract: The mechanical and thermal sensitivity of pyrotechnic flash compositions consisting of mixtures 
of potassium nitrate (KNO3), sulphur (S) and aluminum (Al) with varying particle sizes of KNO; and Al 
indicated that, irrespective of the composition of the cracker mixture, all the compositions were found to 
be thermally and mechanically sensitive. Although the impact sensitivity results reflected the change in the 
surface area of the particle sizes, the changes were within a narrow range of limiting impact energy (LIE) 
(7.5—9.1 J). Further it was difficult to pinpoint a particular sieve fraction as sensitive since the response to 
explosion depended normally on the flash composition and the particle size but also on the density and the 
compactness of the chemicals 


DSC studies on the effect of the Al particle size showed that a decrease in the Al particle size led to a second 
exothermic activity. This behavior should be viewed with caution when considering safety aspects. 


Keywords: pyrotechnics, flash composition, fireworks, impact sensitiveness, friction sensitiveness, 


differential scanning calorimeter, potassium nitrate, sulphur, aluminum 


Introduction 


In recent years frequent incidents have been 
reported involving fireworks during _ their 
processing, storage and transportation.’ Large 
quantities of different categories or types of 
fireworks are manufactured and demand for them 
is rising continuously. In general, the composition 
of fireworks is basically a mixture of sulphur, 
phosphorus, chlorates, nitrates and pure aluminum 
metal powder. These mixtures are highly sensitive 
to temperature, mechanical impact, pressure and 
friction. Knowledge of the thermal stability, auto- 
ignition temperature, impact sensitivity, frictional 
sensitivity and electrostatic sensitivity of these 
mixtures is required to assess their hazard potential 
as well as to make suitable plans for safety during 
processing, storage, and transportation.”* 


Reported information on the sensitivity to thermal, 
mechanical and electrostatic hazards is scarce. 
Recently efforts have been made by Sivapirakasam 
et al. to study the thermal hazards and the impact 
sensitivity of various compositions of cracker 
mixtures.*° The parameters that influence the 
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mechanical and thermal sensitivity are particle 
size, purity, and moisture content of the chemicals 
and ambient conditions. Though there are a few 
reported studies on the effect of particle size on 
the mechanical sensitivity and thermal stability of 
some pyrotechnic mixtures,’ pyrotechnic flash 
composition mixture consisting of KNO;, S and 
Al has not yet been studied. This work attempts 
to report the effects of the particle sizes of KNO, 
(oxidizer) and Al (fuel) in pyrotechnic flash 
compositions on the mechanical sensitivity and 
thermal stability. 


Particle sizes of the ingredients and 
their effects in pyrotechnics 


Pyrotechnic compositions are sensitive to thermal 
and mechanical stimuli. The degree of the hazard 
depends upon the rate of availability of these 
energies in addition to the particle size of the 
components of the mixture, confinement and 
momentum (the force during impact). Each one of 
the factors has a significant role in enhancing the 
overall hazard potential of the mixture. 
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It has been reported® that the ease of ignition is 
greatly dependent on the particle size and surface 
area of the ingredients. Smaller particle sizes create 
large interfacial areas, and increased numbers of 
atoms at the particle interface which on ignition 
lead to a higher heat of reaction. These effects 
however are not actually consistent due to the 
results of various physical and chemical processes 
occurring concomitantly and in many instances 
competitively. This is particularly true where an 
endothermic phase change absorbs some of the 
heat of combustion.” 


When the particles in a pyrotechnic mixture are 
small and jagged, less energy input (via impact, 
friction and thermal stimuli) is needed to produce 
hot spots. This is because energy is localized at 
the stress points. Also, dislocations, cracks and 
other discontinuities in the crystal structure are 
conducive to the formation of hot spots.'° Ignition 
is achieved when the energy released from the 
initiating hot spots is sufficient in quantity and 
is efficiently captured by the adjacent layers or 
grains of blended chemical compositions. '! 


Experimental 
Materials 


The flash compositions were prepared from 
materials purchased from SD’s fine chemicals, 
Mumbai. KNO; and S were of high purity 
(AnalaR) grade and Al was of LR grade. 


Preparation of sieve fractions of KNO, and Al 


The sieve fractions of KNO; and Al were prepared 
using the standard sieves obtained from Jayant 
Scientific Industries, Mumbai, India. Particles of 
KNO, passed through a 63 um sieve mesh and 
collected in a 53 um sieve mesh were termed as 
+53 um sieved fractions. Similarly, five other 
sieve fractions of +63 um, +75 um, +90 um, 
+105 pm, and +150 um were prepared. The details 
of the sieved fractions separated for the study are 
summarized in Table 1. Similarly Al was separated 
into 5 sieved fractions and the details are given in 
Table 2. 


In order to examine the effect of varying the 
particle size of KNO; on the sensitiveness of a 
flash composition, the sieved fraction of KNO; 
was mixed with the other two components i.e., S, 
Al. The particles of S and Al were passed through 
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Table 1 Particle sizes of KNO3. 


Sample No. Particle Size of KNO; 
1 + 53 um (53 to 63 um) 
2 +63 um (63 to 75 um) 
3 +75 um (75 to 90 um) 
a +90 um (90 to 105 um) 
5 +105 um (105 to 150 pm) 
6 +150 um (150 to 200 um) 


Table 2 Particle sizes of Al. 


Sample No. Particle Size of Al 
1 +37 um (37 to 45 um) 
2 +45 um (45 to 53 wm) 
3 +53 um (53 to 63 wm) 
4 +63 um (63 to 90 um) 
5 +90 um (90 to 105 um) 


a 100 mesh brass sieve. Similarly, in order to 
examine the effect of varying particle sizes of Al 
on the sensitiveness of a flash composition, the 
sieved fraction of Al was mixed with the other 
two components i.e. KNO; and S. The particles 
of KNO, and S were passed through a 100 mesh 
brass sieve. Tables 3, 4, 5 and 6 show the various 
flash compositions and the corresponding mixture 
fractions taken for studies. 


Measurement of impact sensitivity 


Impact sensitivity measurements on flash 
compositions with varying KNO; and Al 
particle sizes were carried out according to the 
procedure outlined in the United Nations (UN) 
Recommendations on the transport of dangerous 
goods.’? The design and principle of the equipment 
were similar to those of a drop fall hammer of 
BAM standards. The details of the equipment 
employed in this study are presented elsewhere.’ 
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Measurement of friction sensitivity 


Friction sensitivity measurements on flash 
compositions with varying KNO, and Al particle 
sizes were carried out by BAM (friction tester) 
according to the procedure outlined in the United 
Nations (UN) Recommendations on the transport 
of dangerous goods." 


Thermal studies 


A Differential Scanning Calorimetry (DSC) 
module 2910 model from TA Instruments was 
used for thermal stability measurements for the 
flash compositions with different KNO; and Al 
particle sizes. The studies were conducted with 
a sample size of 2 mg under pure nitrogen gas at 
a flow rate of 100 pl min '. The equipment and 
the experimental conditions employed have been 
reported in detail elsewhere.” 


Results and Discussion 


Effect of KNO; particle size on impact 
sensitivity 


Impact sensitivity testing results for 3 different 
flash composition mixture ratios consisting of 
KNO,, S and Al with varying KNO, particle sizes 
in the range of 53-150 um are shown in Figure 1. 
The plot shows that all the mixture compositions 
studied were impact sensitive and Limiting Impact 
Energy (LIE) was in the range of 7.5 to 9.1 J. 
Hence, these mixtures could be groupedas category 
III explosives according to the classification of 
Andreiev-Beliave.'*'° From Figure 1 it can be 
seen that an increase in the particle size of KNO, 


9.5 
1 -KNO3=65%, S=20%, Al=15% 
2 - KNO3=50%, S=11%, Al=39% 
5 94 3 - KNO3=53%, S=17%, Al=30% 
> 
> 
2 
w 8.54 
2) 
o 
Q 
E 
» 84 
= 
E 
3 
1.554 
7 T T T T 
45 70 95 120 145 


Particle Size of KNO3, mm 


Figure 1 Effect of KNO; particle size of flash 
compositions on the impact sensitivity. 
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decreased the sensitivity to impact initially up to 
63 um and then increased it. This trend was due to 
the fact that impact sensitivity depends not only on 
the flash composition and the particle size but also 
on the particle shape, density and compactness of 
the chemicals. 


Effect of Al particle size on impact sensitivity 


Figure 2 shows the impact on the flash composition 
mixtures with varying Al particle sizes. The flash 
composition mixtures with the lowest particle 
sizes showed high sensitivity compared with 
those with higher particle sizes. This trend could 
be attributed to the increase in surface area as 
the Al particle size decreased. All the mixtures 
studied were impact sensitive and varied within a 
narrow range of LIE (8.5—11.5 J) and hence the 
flash composition mixtures could be grouped as 
category III explosives as per the classification of 
Andreiev-Beliave.'*!° 


Effect of KNO, particle size on friction 
sensitivity 


Figure 3 shows the effect of KNO; particle size on 
the friction sensitivity of various flash composition 
mixtures consisting of KNO3, S and Al. The friction 
sensitivity decreased with increasing KNO, 
particle sizes. The varying particle size of KNO; 
set the lowest friction load of 144 N for the flash 
composition studied. KNO; being an oxidizer, the 
effect of particle size was appreciable. 


1 - KNO,=65%, S=20%, Al=15% 
13 2 -KNO;=53%, S=17%, A=30% 
3 - KNO;=50%, S=11%, AI=39% 


114 


Limiting Impact Energy, J 


35 55 75 95 
Particle Size of Al, um 


Figure 2 Effect of Al particle size of flash 
compositions on the impact sensitivity. 
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Figure 3 Effect of KNO; particle size of flash 
compositions on the friction sensitivity. 


Effect of particle size of Al on friction 
sensitivity 


The effect of Al particle size on the friction 
sensitivity of flash composition is presented in 
Figure 4. Although these mixtures were shown to 
be sensitive to friction there was no appreciable 
change with the change in Al particle size. 
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Figure 4 Effect of Al particle size of flash 
compositions on the friction sensitivity. 


Effect of KNO, particle size on thermal 
stability 


Table 5 shows the influence of KNO; particle size 
on the thermal decomposition of a pyrotechnic 
mixture. This is followed using DSC plots stacked 
according to the varying particle sizes (Figures 5, 
6 and 7). In general, for all the flash compositions, 
the exothermic peak recorded between 299 °C 


Table 3 Effect of KNO; particle size for the various flash composition mixtures on mechanical sensitivity. 


Composition (wt%) Particle Size (tum) 


53 
KNO; = 50% a 
S) = 11% 90 
Al = 39% 

105 

150 

53 
KNO, = 53% ie 
Ss = 17% . 
Al = 30% 

105 

150 

53 
KNO; = 65% _ 
Ss = 20% is 
Al = 15% 

105 

150 
Page 42 


Limiting Impact Energy (J) 


Frictional Limiting Load (N) 


8.4 216 
9.0 240 
9.0 240 
8.8 288 
8.4 360 
8.0 360 
7.7 168 
8.8 144 
8.6 144 
8.5 160 
8.6 252 
8.2 288 
8.6 252 
9.3 324 
9.0 324 
8.9 324 
8.8 360 
8.6 360 
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Table 4 Effect of Al particle size for the various flash composition mixtures on mechanical sensitivity. 


Composition (wt%) 


KNO, 
S 
Al 


50% 
11% 
39% 


53% 
17% 
30% 


65% 
20% 
15% 


Particle Size (tum) 


31 
45 
53 
63 
90 
37 
45 
33 
63 
90 
37 
45 
53 
63 
90 


Limiting Impact Energy (J) 


8.6 
9.0 
9.0 
9.4 
9.4 
9.4 
9.6 
9.8 
9.8 
10 
11.3 
11.3 
11.3 
11.3 
11.9 


Frictional Limiting Load (N) 


252 
260 
288 
288 
288 
288 
360 
360 
360 
324 
360 
360 
360 
360 
360 


Table 5 Effect of KNO; particle size for the various flash composition mixtures on thermal sensitivity. 


Composition (wt%) 


Particle Size 


of KNO; (ttm) 


53 
63 
75 
90 
105 
150 
53 
63 
75 


90 


105 
150 
53 
63 
75 
90 
105 
150 


Temperature 


Onset 
CC) 
300.16 
319.66 
306.78 
227.67 
312.88 
310.47 
224.08 
312.2 
310.38 
313.79 
*331.84 
315.53 
314.63 
314.38 
314.98 
313.03 
314.3 
315.42 
314.27 


Peak Maximum 
Temperature 

CC) 
316.57 
316.57 
320.84 
299.03. 
324.15 
320.36 
315.56 
326.12 
324.74 
328.80 

*333.04 
326.22 
325.59 
327.20 
323.78 
326.82 
326.63 
322.27 
328.63 


* represents the DSC data for the second exothermic activity 


Journal of Pyrotechnics, Issue 23, Summer 2006 


Heat of 
Reaction 


Ig!) 
30.25 
48.93 
85.39 

697.8 
73.51 
41.44 

778.9 

290.6 

116.2 

232.6 

*28.16 
78.03 
71.02 

336.7 

101.2 

223.2 

146.2 
56.65 

195.3 


End Temperature 


320.21 
329.90 
324.06 
314.31 
327.07 
324.28 
318.56 
330.52 
327.06 
331.13 
*338.31 
329.57 
328.48 
330.64 
326.48 
329.8 
329.65 
330.42 
330.92 
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Figure 5 DSC Plots showing the influence of KNO; particle size on thermal decomposition of a flash 


composition (KNO; = 50%, S = 11%, Al = 39%). 


and 333 °C was due to the prominent exothermic 
decomposition process. The concomitant endo- 
thermic peak with a peak maximum of 331 °C 
was due to melting (phase change) of unreacted 
KNO,. Though the degree of exothermic activity 
was different with varying particle size, there was 
no specific trend or influence of particle size on 
the thermal decomposition. 


A close examination of the DSC plot revealed 
that the decomposition process was a result of 
the physical and chemical processes occurring 
concomitantly. It appeared that the initial reaction 
process occurred in the solid phase (220 °C to 
330 °C), until the endothermic melting of KNO;, 
and hence the decomposition reaction could not 
proceed further at this temperature. The possibility 
of the existence of sulphur was remote and hence 
the decomposition process ended at this point. 


The DSC plots (Figure 5) for the flash compositions 
consisting of 50% KNO3, 11% S and 39% Al with 
varying particle sizes show a number of features 
in common, except for the one carried out with 
KNO; of 90 um particle size. The 90 um particle 
size mixture showed an early onset of thermal 
decomposition (227°C) and the heat released 
was quite high as compared to the DSC plots 
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of the remaining mixtures of varying particle 
sizes. Perhaps this oddity is because the flash 
compositions depend also on the particle shape 
and compactness of the chemical components. 
Though the changes in particle size in the mixture 
composition had a definite influence on the thermal 
decomposition, it was difficult to conclude which 
particle size and composition was the most ideal 
for flash composition manufacturing. 


The influence of KNO; particle size on the thermal 
decomposition of a pyrotechnic mixture consisting 
of 53% KNO;, 17% S and 30% AI was studied 
using DSC and the resulting plots were stacked 
according to the varying particle sizes (Figure 6). 
In case of the 53-63 um sieve fraction, there was 
an earlier onset temperature for decomposition 
with the release of a high heat of reaction. 


The influence of KNO; particle size on the 
thermal decomposition of a pyrotechnic mixture 
consisting of 65% KNO3, 20% S and 15% Al was 
also studied using DSC and the resulting plots 
were stacked according to increasing particle size 
(Figure 7). This showed similar features to what 
was seen in the earlier compositions. Though a 
change was observed on increased KNO; particle 
sizes, the effect was not very significant. 
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Figure 6 DSC plots showing the influence of KNO; particle size on thermal decomposition of a flash 


composition (KNO; = 53%, S = 17%, Al = 30%). 


Effect of Al particle size on the thermal 
stability 


Table 6 shows the influence of Al particle size 
on the thermal decomposition of the pyrotechnic 
mixture. The resulting plots were stacked 
according to the varying particle sizes (Figures 8, 
9, 10). The onset of decomposition was generally 


Heat Flow (mv) 


Temperature (°C) 


found to be around 320°C for all particle sizes 
and the first process ended around 330 °C. Here, 
in most cases, a second exotherm was observed 
immediately next to the endotherm at 331 °C. 
The endothermic process recorded at the peak 
maximum temperature of 331 °C was the melting 
of KNO, [reference Merck index]. DSC plots 


200 
Universal V3.7A TA instromonts 


wo 


Figure 7 DSC plots showing the influence of KNO; particle size on thermal decomposition of a flash 


composition (KNO;= 65%, S = 20%, Al = 15%). 


Journal of Pyrotechnics, Issue 23, Summer 2006 


Page 45 


Table 6 Effect of particle size of Al for the various mixtures of flash compositions on thermal sensitivity. 


Cc ition (wt%) Particle Size T the 
omposition (wt% of Al (um) emperature 
37 220.03 
45 319.26 
KNO; =50% 
319.82 
S =11% 53 
*331.54 
Al = 39% 
63 320.62 
90 315.9 
318.71 
a7 
*332.25 
312.10 
45 
KNO, =53% *332.15 
N) =17% a 316.78 
Al =30% *317.28 
63 319.91 
321.54 
90 
*331.67 
37 317.69 
45 311.23 
315.68 
KNO, =65% 53 
*332.21 
N) =20% 
324.6 
Al =15% 63 
*332.66 
325 
90 
*334.95 


Heat of 


Peak Maximum re End Temperature 
Temperature (°C) Ge") (°C) 
332.75 263.2 331.8 
325.68 42.24 329.93 
326.65 27.84 330.77 
*332.28 *11.03 *338.74 
325.84 47.65 329.91 
322.95 34.15 329.65 
326.82 76.7 330.95 
*333.45 *125 *337.49 
329.77 105.5 330.87 
*333.04 *81.43 *345.84 
328.84 101.5 330.82 
"333.13 *23.43 *338.64 
326.24 66.5 330.9 
326.07 37.46 330.81 
*332.68 *17.75 *339.99 
326.35 79.87 329.91 
332.75 1007 353.77 
325.76 88.6 330.91 
*334.29 *103.9 *353.54 
330.41 44.88 331.71 
*338.45 *240.6 *363.17 
331.48 4.84 332.25 
*343.07 *81.81 *353.72 


* represents the DSC data for the second exothermic activity 


also revealed that the decomposition process was 
the result of physical and chemical processes 
occurring concomitantly. The effect of particle 
size seemed to be inconsistent. It appeared that the 
initial reaction was taking place in the solid phase 
(220-330 °C), until the endothermic melting of 
KNO,. The concomitant reaction after KNO, 
melting appears to proceed in the liquid phase. 


Figure 8 shows the DSC plot for flash compositions 
consisting of 50% KNO3, 11% S and 39% AI with 
varying Al particle sizes. The second concomitant 
reaction could not be observed, because the first 
reaction could not provide enough enthalpy for 
the second reaction to proceed due to absorbance 
of the heat by the endothermic phase change. 
However, in some cases the second exothermic 
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reaction could proceed. From Table 6, it was 
observed that the onset temperature increased with 
increase in particle size and remained constant. 
The decrease in the Al particle size increased the 
heat of reaction. 


The influence of Al particle size on the thermal 
decomposition of pyrotechnic mixture consisting 
of 53% KNO3, 17% S and 30% AI was studied 
using DSC and the resulting plots were stacked 
according to increasing particle size (Figure 9). 
In all the DSC runs with different particle sizes, 
the endothermic melting of potassium nitrate with 
concomitant second exothermic activity followed 
the first prominent exothermic peak. 


The influence of Al particle size on the thermal 
decomposition of pyrotechnic mixture consisting 
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Figure 8 DSC plots showing the influence of Al particle size on thermal decomposition of a flash 
composition (KNO; = 50%, S = 11%, Al = 39%). 
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Figure 9 DSC plots showing the influence of Al particle size on thermal decomposition of a flash 
composition (KNO; = 53%, S = 17%, Al = 30%). 
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Figure 10 DSC plots showing the influence of Al particle size on thermal decomposition of a flash 


composition (KNO; = 65%, S = 20%, Al = 15%). 


of 65% KNO;, 20% S and 15% Al was studied 
using DSC and the resulting plots were stacked 
according to increasing particle size (Figure 10). 
The trend of the curve was similar to the two 
previous compositions except for the one with 
particle size of 45 um. It was very interesting 
to note that the concentration of KNO, was the 
highest in the mixture and that all of them were 
involved in the thermal decomposition processes. 
As KNO, was the oxidizer in the flash composition 
a maximum amount of heat release was expected. 
The results are complex and it is difficult to 
interpret the role of each component in promoting 
the decomposition reaction. 


Effect of Particle Size of KNO, 
(oxidizer) and Al (fuel) on Process 
Safety of Flash Composition 


The mechanical sensitivity analysis showed that the 
decrease in particle size of KNO; had an adverse 
effect on sensitivity (both impact and friction) 
to a greater extent than the Al particle size had. 
However, the flash compositions studied could 
be categorized as class III explosives sensitive to 
impact. 
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The DSC studies showed that the decrease in 
particle size of KNO, and Al set the lowest 
onsets at around 220°C for the first exothermic 
activity. The chance of thermal hazard during the 
processing of flash compositions below the onset 
temperature of DSC, though remote, had to be 
so confirmed through the adiabatic calorimetric 
test. Given a margin of 100 °C, it would be stated 
that below 120 °C, the flash composition mixture 
would be stable provided there was no source of 
ignition or mechanical (impact/friction) effects of 
the flash composition. 


The DSC studies on the effect of Al particle size 
showed that the decrease in Al particle size leads 
to a second exothermic activity. This behavior 
needs to be viewed with caution from the point 
of view of safety. The heat content from the first 
exothermic activity itself is enough to undergo 
good flash reaction by the mixture,° whereas the 
heat release due to the second exothermic activity 
is undesirable for the flash composition as this may 
lead to several cascading explosions. However, in 
order to improve the quality of cracking of flash 
composition, firework manufacturers generally 
tend to use finely divided Al powder in greater 
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quantities than necessary. DSC studies revealed that 
the use of finer Al powders in a flash composition 
would only lead to secondary exothermic activity 
culminating in accident situations of higher 
magnitude. This secondary exothermic activity 
can lead to secondary explosions in any hazardous 
situation during processing and storage of flash 
compositions as the quantities are going to be 
high. 


Conclusions 


Mechanical and thermal studies on the sensitivity 
of pyrotechnic flash compositions with varying 
KNO, and AI particle sizes indicated that all 
the flash compositions studied were found to be 
sensitive. As the KNO, particle size increased, 
the sensitiveness to impact initially decreased up 
to 63 um and then increased. This was due to the 
fact that impact sensitivity depended not only on 
the flash composition and particle size but also on 
the particle shape, density and compactness of the 
chemicals. The friction sensitivity decreased with 
the increase in KNO; particle size. The lowest 
Al particle size exhibited high impact sensitivity 
compared with the higher particle sizes. But there 
was no appreciable effect on friction sensitivity 
with change in Al particle size. 


The DSC studies on the effect of Al particle size 
showed that the decrease in Al particle size led to a 
second exothermic activity. Hence, the use of finer 
Al powders in a flash composition would only lead 
to secondary exothermic activity. This secondary 
exothermic activity can lead to secondary 
explosions in any hazardous situation during 
processing and storage of flash compositions. 
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Interrelation Between Impact, Friction And Thermal Energy 


In A Pyrotechnic Flash Reaction 
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Abstract: Firework chemical mixtures are sensitive to thermal and mechanical stimuli and lead to 
many explosive incidents. Experimentally determined thermal and mechanical (impact and friction) 
sensitivity data of a flash composition mixture are subjected to statistical and graphical analysis in 
order to understand the mechanism of triggering accidents. The interrelationship study reveals that 
irrespective of the nature of stimuli, explosion is the final event and occurs due to a thermal mechanism. 
This study shows that under severe impact thermal stimuli can occur. If the thermal stimuli are equal to 
or greater than the activation energy of the composition then ignition of the flash composition will occur. 


Keywords: flash composition, fireworks, mechanical and thermal sensitiveness, correlation analysis 


Introduction 


During the firework manufacturing process, 
chemicals are initially mixed to produce a 
reasonably homogeneous blend of oxidizer, fuel, 
colour enhancing chemicals and binders. During 
these operations impact, friction, spark and heat 
stimuli may occur, and under certain conditions 
one or more of these stimuli may be enough to 
cause ignition of the composition. The sensitivity 
of the flash composition to these stimuli depends 
upon the chemical components, purity, particle 
size, moisture content and packing density. It is 
therefore extremely important to understand the 
sensitivity of the chemical mixture to external 
stimuli. 


Studies on thermal stability,' impact sensitivity 
and friction sensitivity” of firework compositions 
have been reported. The reported information 
cannot be used directly for determining the safety 
limits for storage, processing, and transportation 
of firework compositions because no attempts 
have been made to study the thermal, mechanical 
and electrostatic hazards together to pinpoint the 
reasons for accidents to occur. However, in a 
few reported studies,*© attempts to correlate the 
mechanical initiation of organic high explosives to 
kinetics of thermal decompositions are evident. 


Ho and Fong* compared the impact sensitivity 
of the various propellants with the thermal 
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decomposition data at 20°C min and showed 


that the impact energy had a good correlation 
with thermal decomposition for the composite 
propellants with the same binder : oxidizer weight 
ratio. The results at 20°C min! were used for 
comparison so as to emphasize the thermal effect 
most pronounced under rapid heating conditions. 


Wenograd* suggested that impact sensitivities of 
organic high explosives were governed by the 
thermal decomposition processes which took 
place at the widely varying temperature generated 
under the impact mass. He found that the 
temperature at which an explosion would occur 
within 250 microseconds (a time comparable to 
the interval under the impact mass) varied greatly 
among explosives. 


Bowden> made efforts to study the mechanism 
of impact initiation. The authors suggested that 
initiations stemmed from hot spots in the explosive 
mass generated by a number of possible routes 
including viscous heating, frictional heating and 
adiabatic compression of entrapped gases. They 
concluded that to cause fires in PETN and NG, 
these hot spots must reach temperatures of at least 
430-500 °C. 


Field et al.° suggested shear banding as a possible 
source of hot spot formation on impact. In this 
mechanism, localized plastic flow of material 


Page 51 


following its structural collapse under compressive 
stress was related to hot spot formation. 


Of all the firework compositions, flash 
compositions are real explosives that detonate, 
if a sufficient quantity of powder is present in 
bulk form, even if unconfined.’ It is necessary to 
study the thermal, impact and friction sensitivity 
of these compositions so that the interrelation and 
the functioning of mechanical initiation to thermal 
decomposition can be understood. The work 
reported here is focused on this objective. Further 
in this paper, for the first time, Pearson’s correlation 
analysis’ has been employed to interrelate the 
mechanical and thermal decomposition parameters 
and the resulting data have been analyzed using 
graphical methods. 


Experimental 
Materials 


The chemicals used for the preparation of the 
flash composition were of commercial grade and 
obtained from a firework manufacturing company 
situated in the southern part of Tamilnadu, India. 
The purity and assay of the chemicals were: KNO3: 
91.6%, S: 99.84% and Al: 99.71%. The chemicals 
were passed through a 100 mesh brass sieve. The 
samples were stored in an airtight container and 
kept away from light and moisture. 


The mixture compositions varied in the range of 
50-65% potassium nitrate, 5—20% sulphur, and 
45-15% aluminum. The range was kept wide so 
as to cover the range of compositions employed 
in different fireworks industries in and around 
Tamilnadu, India. 


Thermal studies under isothermal conditions 


Intensive studies with the pyrotechnic flash 
compositions consisting of KNO;, S and Al 
were carried out using Differential Scanning 
Calorimetry (DSC) to analyze their thermal 
stability and understand the importance of the 
role of varying proportions and particle size 
of pyrotechnic flash compositions in inducing 
cracking characteristics. 


DSC module 2910 model (TA Instruments) was 
used for measurement of thermal stability for the 
different flash compositions. The studies were 
conducted with a sample size of 2 mg under 
pure nitrogen gas at a flow rate of 100 1 min '. 
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The equipment and the experimental conditions 
employed have been reported elsewhere. ' 


Thermal studies under adiabatic conditions 


Thermal study of the flash compositions under 
adiabatic conditions was carried out using an 
Accelerated Rate Calorimeter (ARC). An ARC 
1000 supplied by CSI of Austin, TX was used.® 
About | g of sample was loaded into a titanium 
bomb calorimeter and its temperature raised 
incrementally by 5°C min! under heat—wait— 
search mode, until a measurable rate of exothermic 
activity was detected (0.02 °C min‘) or the final 
temperature was attained without any positive 
thermal input. 


Measurement of impact sensitivity 


The impact sensitivity measurements of the flash 
compositions were carried out according to the 
procedure outlined in the United Nations (UN) 
Recommendations on the transport of dangerous 
goods.° The design and principle of the equipment 
were similar to those of a drop fall hammer of 
BAM standards. The details of the equipment 
employed have been presented elsewhere.” 


Measurement of friction sensitivity 


The friction sensitivity measurements of the flash 
compositions were carried out by BAM (friction 
tester) according to the procedure outlined in the 
United Nations (UN) Recommendations on the 
transport of dangerous goods." 


Correlation analysis’ 


If two variables vary such that change in one 
variable affects the change in the other variable, 
the variables are correlated. The degree of 
relationship between the variables under 
consideration is measured through correlation 
analysis. The measure of correlation is called the 
correlation coefficient or correlation index. Thus, 
correlation analysis refers to the techniques used 
in measuring the closeness of the relationship 
between the variables. Although there are several 
methods"! of analyzing the correlations of physical 
and chemical sensitiveness, Pearson’s coefficient 
of correlation is simple and highly reliable for 
measuring the degree of relationship between two 
variables. The correlation coefficient p between 
two random variables X and Y is given as follows 
(equation 1): 
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_ nXXY —LXEY 
(n=X? —(ZX)’)(nzY* -(ZY)*)'”? (1) 


fe) 


The value of the correlation coefficient p always 
lies between +1 and —1. A value of p = 0 indicates 
no correlation. If the value of p is near +1 then 
the variables X and Y are said to be positively 
correlated and if the value of p is near to —1 then 
the variables X and Y are said to be negatively 
correlated. 


Results 


Flash compositions under isothermal 
conditions 


The results of the experiments conducted using 
DSC for the different flash compositions are 
presented in Table 1. The DSC plots of the flash 
compositions for varying sulphur concentration 
are shown in Figure 1. It can be seen that with 
increasing sulphur content, the onset temperature 


Table 1 Explosive parameters for the various flash compositions. 


q 0 Peak Friction 
oe Mixture Components (wt7%) ta (C) Temperature AH(Jg') LIE(J) © sensitivity 
KNO, S Al CC) N 
1 65 20 15 310.37 319.69 246.2 6.5 360 
2 65 17 18 310.68 314.79 253.2 7.9 240 
3 65 14 21 313.98 321.79 466.7 7 240 
4 65 5 30 305.85 309.37 53.83 7.9 360 
5 62 20 18 312.76 321.09 332.1 6 216 
6 62 17 21 311.83 320.92 252.5 7.8 216 
7 62 14 24 309.85 314.13 56.27 7.2 240 
8 62 11 27 308.14 314.02 79.78 8.8 192 
9 60 75 32.5 304.06 308.14 33.94 6.5 216 
10 59 20 21 310.22 314.49 36.24 7.6 216 
11 59 17 24 310.5 319.5 358.8 7.6 180 
12 59 14 27 303.58 321 966.2 6.1 192 
13 56 20 24 312.91 320.5 525.8 6 216 
14 56 14 30 306.35 320.3 397.2 6.2 192 
15 53 20 27 312.11 313.56 372.3 7.0 216 
16 53 17 30 310.37 318.24 409.9 6.3 216 
17 53 14 33 311.84 313.88 346 6.4 240 
18 53 11 36 310.59 312.87 227.6 6.6 240 
19 52.5 75 40 307.5 311.2 53.64 7.2 240 
20 50 20 30 308.61 319.81 335.6 6.7 216 
21 50 17 33 310.97 321.31 534.1 5.9 288 
22 50 12.5 37.5 306.75 318.82 367.9 6.0 288 
23 50 11 39 311.61 315.3 326.3 6.1 360 
Table 2 Summary of ARC data of flash composition. 
There Onset Final Adiabatic _ Absolute — Heat of 
Inertia (0) Temperature Ty Temperature T; temperature rise temperature rise reaction AH, 
CC) (°C) ATC) AT CC) (jg ) 
4.84 191 450 259 1253.6 1311.25 
Journal of Pyrotechnics, Issue 23, Summer 2006 Page 53 


2.0 = ees 
1-S-11%, Al-27% 
2 - S-17%, Al-21% 


3-S -20%, Al-18% . 
321.09 C 


1.5 
320.92°C 


= 
oO 


313.95°C 


= 
a 


Heat Flow, W/g 


0.0 


é 
| eSieeeees eee ee 
oe Cee e.nee wns 


2 aes ‘ed 
ot 
WRT Oe ene sews on we oe teen cm! <= <"= 


220 240 260 280 300 320 340 360 380 


Universal V2.7A TA Instruments 
° 
Temperature, C 


Figure 1 DSC plots of pyrotechnic flash composition of varying sulphur content (KNO;: 62% fixed). 
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Figure 2 Se/f-heat rate vs. temperature plot for thermal decomposition of flash composition (KNO;: 8: 
AL 53-17: 30). 
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for exothermic decomposition advances to a higher 
value to increase the heat of decomposition. 


Flash compositions under adiabatic conditions 


The self-heat rate plot for thermal decomposition 
of flash composition consisting of potassium 
nitrate, sulfur, aluminum in the ratio of 53 : 17 : 30 
is shown in Figure 2. ARC data are summarized 
in Table 2. 


Mechanical sensitiveness of flash compositions 


The results of mechanical sensitiveness (impact 
and friction) are shown in Table |. The limiting 
Impact Energy (LIE) falls in the range of 5-8 J 
for the compositions studied, which may be 


categorized as class III explosives according to the 
1965 classification by Andreieve-Beliaev.” The 
friction-limiting load falls in the range of 192— 
240 N for the compositions studied. The impact 
energy and friction-limiting load vary when the 
concentration of any one of the components of the 
mixture is changed. This behavior is due to the 
sensitivity and reactivity of each component. 


Correlation analysis on mechanical 
sensitiveness and thermal decomposition 


The results of thermal, impact and_ friction 
sensitivity data (Table 1) were subjected to Karl 
Pearson’s correlation analysis to understand the 


Table 3 Correlation of mechanical sensitiveness and thermal decomposition parameters. 


Variables 


Limiting impact energy vs. onset temperature 
Limiting impact energy vs. peak temperature 
Limiting impact energy vs. heat of reaction 
Friction sensitivity vs. onset temperature 
Friction sensitivity vs. peak temperature 
Friction sensitivity vs. heat of reaction 


600 


500 


400 


300 


Heat of Reaction, J/g 


200 


100 


11 13 15 


Correlation coefficient 


Significance 


+0.26 A weak positive correlation 
—0.6 A strong negative correlation 
—0.6 A strong negative correlation 
+0.23 A weak positive correlation 
—0.13 A weak negative correlation 
+0.17 A weak positive correlation 


1 - KNO3= 50% 
2 - KNO3= 53% 


17 19 21 


Wt. % of Sulphur 


Figure 3 Effect of sulphur concentration on the heat of decomposition of flash compositions. 
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interrelation between the mechanical and thermal 
explosive sensitiveness. 


Equation (1) was employed to determine the 
correlation coefficient p between impact friction 
and thermal sensitivity. The results are summarized 
in Table 3. 


Discussion 


The plot between the decomposition energy 
and sulphur concentration (Figure 3) showed 
that with increasing sulphur concentration, the 
decomposition energy (H decomposition) release 
increased and reached a maximum value at 17% 
concentration. Above this, the H decomposition 
decreased, perhaps due to an increase in 
concentration of other two components. Table 1 
show that the flash composition was found to 
have explosive characteristics between 11 and 
20% of sulphur concentration. Thus, sulphur 
concentration in the flash composition appeared 
not only critical, but should also be around the 
optimum level to exhibit good flash properties. 
The onset temperature in DSC was above 303 °C 
(Table 1). However, the reactive potential of 
flash composition under adiabatic conditions was 
severe, and a lower onset for decomposition was 
recorded as 191 °C (Table 2). A peak self-heating 
rate of 2.625 °C min ' was registered at 302 °C. 
The adiabatic temperature rise for the process 


Heat of Reaction, J/g 


was 260°C. Under adiabatic conditions flash 
compositions decomposed slowly until 250°C 
(1700 min) (Figure 2), and beyond this the rise in 
temperature was sudden and sharp. This showed 
that, under adiabatic conditions, flash composition 
underwent vigorous decomposition. 


Interrelation between limiting impact energy 
and thermal decomposition 


The results of correlation analysis are summarized 
in Table 3. The strong negative correlation between 
LIE and AH and peak temperature is an indication 
of the fact that a pyrotechnic mixture, when 
subjected to an impact force, is triggered to an 
energetic response i.e. it explodes. In Figure 4 the 
interrelation between LIE and Ad is graphically 
shown. It can be seen from the graph that, within 
the experiments conducted, AH increases with the 
concentration of KNO; to a maximum of 59%; 
between 59 and 62% AH drops to its minimum; 
beyond 62% KNO; AH increases again. Moreover, 
with increase in concentration of KNO3, the impact 
energy initially decreases and later there is a rapid 
increase in impact energy before it decreases 
again. A close examination of the graph reveals 
that there is an inverse relationship between LIE 
and AH. That is, at lower impact energy (higher 
sensitivity), AH released is higher. This shows that 
the degree of energetic response is dependent on 


Limiting Impact Energy, J 


59 


Wt. % of Potassium Nitrate 


—a— Heat of Reaction vs Wt. % of Potassium Nitrate 
—#-— Limiting Impact Energy vs Wt. % of Potassium Nitrate 


Figure 4 /nterrelation between heat of reaction and impact energy of flash composition at fixed S = 14 


wt%. 
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325 


320 


315 


310 


Peak Temperature, °C 


305 


300 


53 56 


59 


Limiting Impact Energy, J 


62 


65 


Wt. % of Potassium Nitrate 


—m— Peak Temperature vs Wt. % of Potassium Nitrate 
—¢— Limiting Impact Energy vs Wt. % of Potassium Nitrate 


Figure 5 /nterrelation between peak temperature and limiting impact energy of flash composition at fixed 


S= 14%. 


the concentration of mixture constituents and not 
the impact force. 


In Figure 5, the interrelation between LIE and 
peak temperature of flash compositions at a fixed 
sulphur concentration (14 wt%) is graphically 
shown. It can be seen that higher impact sensitivity 
leads to higher peak temperature. 


From Table 3, the weak positive correlation of 
the impact energy to onset temperature can be 
attributed to the fact that the temperature generated 
through impact energy should be more than that of 
the onset temperature for the ignition to occur. In 
the impact sensitivity experimental measurement, 
the impact energy is measured when the explosion 
occurs. It is noted that the time factor between the 
applied impact force and explosion occurring is 
sudden and almost instantaneous, and may be of 
the order of microseconds (up to 250 microseconds 
for most high explosives‘). In the thermal treatment 
of the composition under adiabatic conditions, it 
is seen that sudden and instantaneous explosion 
is achievable only when the temperature on the 
mixture is raised to a temperature beyond 302 °C 
for the flash compositions (self-heat rate plot 
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Figure 2). 


Therefore, it is hypothesized that in impact 
sensitivity measurement, the temperature 
generated under impact mass should be more 
than the temperature of 302 °C for the sudden and 
instantaneous explosion to occur. However, this 
does not mean that explosion will not occur when 
the impact force produces a temperature lower 
than 302°C (under conditions of less impact). 
Flash composition is known to undergo a sort 
of self-heating type of explosive decomposition 
as observed in the ARC studies. In a practical 
situation, where the impact force can produce an 
initial temperature of 191 °C (development of hot 
spots) the onset of explosion may occur, however 
there is some induction time before it reaches the 
critical temperature of 302°C. Therefore, any 
impact force that contributes to rise in temperature 
of flash composition around 191°C or more 
is certainly dangerous. Further, the vigor and 
induction time of the explosion primarily depend 
on the compactness, density and particle size of 
the mixture. 


There is no direct correlation between thermal and 
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Instantaneous 
Explosion 


Pyrotechnic Flash 
Composition 


HOT SPOT DEVELOPS 


Spark 


Ignition 


If Temperature 
< 302°C 


Exothermic 
Reaction 
Propagation 


Figure 6 Flow chart showing the initiation mechanism for explosive decomposition of pyrotechnic flash 


composition. 


impact sensitiveness, to either predict one from 
the other or to predict which of these forces can 
come together to trigger a thermal explosion. It 
is hypothesized that impact stimuli cause thermal 
stimuli for the flash composition to undergo thermal 
explosion. Under severe impact, thermal stimuli 
can occur immediately and lead to a catastrophic 
thermal explosion. Irrespective of the nature of 
the stimulus, explosion occurs through thermal 
mechanism only. This means that the impact or 
other kind of stimulus can only initiate the thermal 


Page 58 


mechanism by providing the minimum threshold 
energy needed/necessary to a reaction temperature 
of 191 °C observed experimentally as the onset 
point for thermal explosion in ARC. Figure 6 
shows the initiation mechanism for explosive 
decomposition of a pyrotechnic flash composition 
arising from various stimuli. It is thus possible to 
interrelate the mechanical form of energy leading 
to the threshold energy (AE) observed in the ARC. 
This provides a means of suggesting a predictive 
correlation in such explosive systems. The degree 
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of explosivity also depends on other factors such 
as compactness, particle size and shape and other 
environmental conditions. 


The impact energy has generally been considered 
to be insufficient to heat the whole pyrotechnic 
charge to adiabatic temperature. This may mean 
that the initiating mechanism is most likely to 
be manifested under impact pressure due to the 
following factors: 


1. Adiabatic compression of trapped gas. 


2. Viscous heating of material rapidly extruded 
between impacting surfaces or grains. 


3. Friction between impacting surfaces, the 
explosive materials and/or grit particles in the 
explosive layer. 


4. Localized adiabatic shear of the materials 
during mechanical fail. 


Correlation for the flash composition shows 
that there exists a relation between LIE, AH and 
peak temperature. It can be concluded that flash 
composition mixtures sensitive to impact will also 
be thermally sensitive and the higher the impact 
sensitivity, the higher the heat of reaction will be. 
The degree of energetic response will primarily 
depend on the mixture compositions and not on 
impact pressure. A minimum of 302 °C is required 
for the instantaneous thermal decomposition under 
impact pressure for a flash composition. However, 
any impact force which contributes to a rise in the 
temperature of flash composition around 191 °C 
or more is certainly dangerous. This is because 
flash composition mixture is known to undergo a 
self-heating type of explosive decomposition as 
observed in ARC studies. 


Interrelation between friction sensitiveness 
and thermal decomposition 


It can be seen from Table3 that the friction 
sensitivity and thermal decomposition are weakly 
correlated. It is difficult to offer any scientific 
explanation at this stage; an acceptable reason 
may be the lack of precision in the measurement 
of friction sensitiveness. In the BAM friction 
sensitiveness measurement apparatus (employed in 
this study) the pyrotechnic mixtures are subjected 
to a localized frictional load and do not cover the 
entire sample subjected to the test. Considering 
the physical nature of the sample and the expected 
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chemical mechanism available for reaction, the 
friction sensitivity data obtained may not be a true 
representation of real life situations. Further work 
is in progress. 


Conclusions 


The correlation analysis has proved that there 
exists a relation between impact and thermal 
sensitiveness for a pyrotechnic flash composition. 
The inverse relationship between limiting impact 
energy and heat of reaction shows that the 
degree of energetic response is dependent on the 
concentration of mixture constituents and not the 
impact force. Higher impact sensitivity leads to 
higher peak temperature during thermal stimuli. 
The correlation analysis has also predicted that flash 
composition mixture sensitive to impact will also 
be thermally sensitive. Further it is hypothesized 
that impact stimuli cause thermal stimuli for the 
flash composition to undergo thermal explosion. 
Under severe impact, thermal stimuli can occur 
immediately, leading to ignition of compositions. 
Thus through this correlation analysis, a 
satisfactory explanation could be projected for an 
accident triggering mechanism. Further research 
is required to relate the numerical factors (onset 
temperature, heat of reaction, limiting impact 
energy) to gain a thorough understanding of the 
accident triggering mechanism. 
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Abstract: A pyrotechnic flash composition consisting of 53% KNO3, 30% Al and 17% S is subjected to 
Accelerating Rate Calorimetry (ARC) studies. The onset point for thermal explosion is 191 °C resulting 
in the generation of a considerable quantity of gaseous products. The mixture is vulnerable to thermal 
hazards. There is good agreement between the predicted and experimental self-heat rates determined 


using adiabatic thermo kinetics. 


Keywords: Flash composition, ARC, fireworks, adiabatic thermo kinetics, thermal characterization 


Introduction 


In recent years, frequent accidents during 
processing, storage and transportation have been 
reported in the fireworks industry.' This is of 
great concern because large quantities of different 
types of fireworks are manufactured in India and 
demand for them is steadily increasing. Generally, 
the composition of fireworks is a mixture of 
oxidizer, fuel, igniter, binder and color enhancing 
chemicals. These mixtures have high sensitivity 
to temperature, impact, friction and electrostatic 
stimuli. A thorough knowledge of thermal stability, 
auto-ignition temperature, impact sensitivity, 
frictional sensitivity and electrostatic sensitivity 
of these mixtures is imperative to assess the 
hazard potential.” Also, it should lead to a suitable 
plan for safety during processing, storage and 
transportation. 


Chemical reactions of pyrotechnics produce large 
amounts of heat when confined to a closed system 
and result in thermal explosion. Although there 
are numerous thermal measurement techniques 
available to characterize the hazardous nature 
of pyrotechnic mixtures, Accelerating Rate 
Calorimetry (ARC) is the only adiabatic and 
versatile calorimetry that produces reliable data. 
Because ARC measurements are conducted 
adiabatically (7.e. no heat losses), the result can be 
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effectively correlated with the behavior of energetic 
materials in bulk. The information obtained from 
ARC experiments relates to the onset temperature, 
self-heat rates and pressure activation energy for 
an exothermic reaction. The ARC data can be used 
to set ceiling temperatures and pressures for safe 
operation, storage and transportation. 


In the past, researchers have studied the thermal 
stability and kinetics of pyrotechnic mixtures 
using Differential Scanning Calorimetry (DSC).* 
The thermal data obtained from DSC could not be 
used for determining safe operating temperatures 
due to the uncertainties associated with the very 
small quantity of samples (2-5 mg) used in the 
experiments and poor reproducibility ofresults and 
non-adiabatic experimental conditions. Subjecting 
flash composition mixtures to ARC studies would 
throw light on the behavior of these samples 
under adiabatic conditions; i.e., under conditions 
of bulk storage, handling and _ transportation. 
Such a study has not been attempted, except for 
a theoretical paper detailing the suitability of 
ARC for studying the thermal decomposition of 
pyrotechnic mixtures.* In the present study the 
thermal data from ARC and the thermo kinetics of 
a pyrotechnic flash mixture consisting of potassium 
nitrate (KNO;), sulphur (S), and aluminum (A1) 
have been studied. 
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Experimental 


Materials: Preparation of flash composition 
mixture 


The chemicals used in this study were of 
commercial grade and obtained from a fireworks 
chemical manufacturing company situated in 
southern Tamilnadu, India. The purity and assay 
of the chemicals were KNO,: 91.6%, S: 99.84% 
and Al: 99.1%. The flash composition consisting 
of potassium nitrate, sulphur and aluminum in the 
ratio of 53 : 17 : 30 was mixed using a wooden 
spatula in a non-flammable container, and each 
time a sample size of | g was prepared. The sample 
was then stored in an airtight container and kept 
away from light and moisture sources. 


Method: Accelerating Rate Calorimeter 
(ARC) experiments 


The ARC used in this study was an ARC 1000 
supplied by CSI of Austin, TX. The working 
principle, design description, and operational 
details of ARC are well cited in the literature.° 
ARC measurements were made using a titanium 
sample vessel in heat—wait-search mode. Before 
loading of the sample, the bomb was flushed with 
inert nitrogen gas and precautions were taken not 
to allow air to enter during the sample loading as 
well as during attachment of the sample vessel 
to the instrument. After connection, the sample 
vessel was pressurized to 2500 psi nitrogen gas 
to ensure that there was no leak and that the air 
in the assembly was replaced. The instrument was 
switched to step mode at an initial temperature of 
80 °C, and a wait time of 15 min was set prior to 
entering the search mode. About | g of sample was 
loaded into the titanium bomb of the calorimeter, 
and its temperature was raised incrementally 
by 5°C min ' in heat—wait-search mode, until 
a measurable rate of exothermic activity was 
detected (0.02 °C min’') or the final temperature 
was attained without any positive thermal input. 
The self-heat rate, time, temperature, and pressure 
data were obtained as ARC output. 


Overview of adiabatic thermo kinetics® 


The first assumption in the interpretation of 
ARC experimental data is the representation of 
concentration in terms of temperature differences. 
The equivalence of temperature and concentration 
for a simple well-defined chemical reaction is 
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established using the ratio: 
GC _U2F eas (1) 
C. f-7 ar 


where C is the concentration of the reacting 
substance and T is the temperature. The subscript 
0 indicates some initial condition, and F a final 
state in which the substance has been consumed. 
Then AT = T; — Tp is the temperature rise for the 
reaction. It is also equal to the ratio of enthalpy 
to average specific heat. In this relation the 
disappearance of the reacting species produces a 
proportionate increase in the heat energy. The heat 
of reaction, AH can be calculated from 


AH =mCpAT 


where Cp is the average heat capacity, and m is 
the mass of the sample. 


The heat generated in an exothermic reaction is 
used in three ways viz., to heat the material, the 
container or bomb and the surroundings. The heat 
being used up in heating the sample mass depends 
on the specific heat. The proportion of heat used in 
heating the container is called thermal inertia (0), 
which is expressed as ¢ = [heat capacity of sample 
(S) and container or bomb (B)]/[heat capacity of 
sample]. 


_ Ss: ps. 
*-~mCn 
galt mo op 
mC, 
(2) 


Incorporating the effects of thermal inertia (@), the 
corrected heat of reaction AH, is calculated using 
equation (3): 


AH, =omC pAT (3) 


The question that is basic to the study of the 
relationship of time to explosion is the measure- 
ment and extrapolation of data. Extrapolation 
must involve a concept of concentration since no 
material can continue to self-heat forever. The 
time dependence of concentration for an Nth order 
reaction rate is expressed as follows: 
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= aC _ yon (4) 


dt 


where Cis the concentration, kis the rate coefficient 
and ¢ is the time. When equations (1) and (4) are 
used, additional temperature dependence appears. 


dC _¢ (At) =o. 


dt dtl AT AT dt 
a es a ee 
ee CN-LLAT 5 
ee at [ AT " ©) 


Here mr is defined as the rate of temperature 
increase (or slope of the graph of T vs. #), i.e. the 
self-heat rate. To remove this extra temperature 
dependence, a modified rate is defined as the 
pseudo rate constant, k*. It is defined in such a 
way that its dimensions for any order reaction are 
reciprocal of time. 


N 
Rsk S= a) 


AT | T, =T (6) 


In practice, k* is evaluated from experimental data 
using the right hand side of the expression. With 
the proper choice of N, &* has the same temperature 
dependence as k and yields a straight-line graph. 


The Arrhenius relationship for determining the 
rate coefficients is 


k= acer] ) 


where 7 is the absolute temperature in Kelvin, 
E is the activation energy, R is the universal gas 
constant and A is the pre-exponential factor. The 
In k* vs. 1/T plot yields a straight line with the 
proper choice of N. The activation energy and 
the pre-exponential factor are calculated by the 
following expressions: 


1-T ky (8) 

k AE 

In A =In— + — 
60 RT (9) 
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Results and Discussion 
Flash composition under adiabatic conditions 


The self-heat rate plot for thermal decomposition of 
flash composition consisting of potassium nitrate, 
sulfur, aluminum in the ratio of 53:17:30 is 
shown in Figure | and the results are summarized 
in Table 1. The onset for reaction occurred at 
191 °C and extended until 450 °C. A maximum 
self-heating rate of 2.625 °C min' occurred at 
302 °C. The adiabatic temperature rise for the 
process was 259 °C. Under adiabatic conditions 
flash compositions decomposed slowly until 
1700 min (250 °C) (Figure 2) and beyond this 
the rise in temperature was sudden and sharp as 
indicated in Figures | and 2. The reaction process 
was also accompanied by a considerable pressure 
rise (Figure 3); the peak pressure observed was 
34 psi (2.312 bar) at 450 °C. As per equation (3), 
the heat of reaction was calculated as 1.311 x 10° 
J g', which is found to be more than the value of 
409 J g' obtained under isothermal conditions.? 
The ARC data showed that the pyrotechnic 
decomposition process under adiabatic condition 
was vigorous and therefore dangerous. 


Thermo kinetics of flash composition 


First order model (i.e. N= 1) kinetics were assumed 
for the decomposition of pyrotechnic flash 
composition. For NV = 1 equation (6) becomes 


10 
1 
c= 
£ 
fo 
2 
fc 0.1 
& 
® 
<= 
0.01 


150 200 250 300 350 400 450 
Temperature, °C 


Figure 1 Se/f-heat rate vs. temperature plot 
for thermal decomposition of flash composition 
(KNO;: S: Al; 53 : 17 : 30). 
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Table 1 Summary of ARC data of flash composition. 


‘These Onset Final 
Peer Temperature 7, Temperature 7; 
inertia (D) (‘C) ° 

4.84 191 450 
k*=k=m{(Tp— T) (10) 


Pseudo rate constants (k*) were calculated using 
equation (10). Then In k* versus the inverse of 
temperature was plotted and the plot obtained 
is shown Figure 4. The straight line obtained 
confirms the assumption that the flash composition 
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Figure 2 Time vs. temperature plot for thermal 
decomposition of flash composition 
(KNO;. S: Al; 53 : 17: 30). 
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Figure 3 Temperature vs. pressure plot for ther- 
mal decomposition of flash composition 
(KNO;. 8: Al; 53 : 17: 30). 
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Adiabatic Heat of 
temperature rise as hey a oc reaction AH, 
AT(C) i Gg’) 

259 1253.6 1.311 x 10° 


mixture follows first order kinetics. 


The slope of the plot is equal to AE/R. As 
per equation (8), the activation energy was 
calculated as 63.99 kcal mol! (268 kJ mol'). It 
is seen that the activation energy obtained under 
adiabatic conditions is close to those found under 
isothermal conditions (199.7 kJ mol ') reported 
by us elsewhere.* This shows decomposition 
under isothermal and adiabatic condition operates 
on the same mechanism. Using equation (9), 
the pre-exponential factor was evaluated as 
6.13 x 10°. Thus the Arrhenius rate law for 
thermal decomposition of flash composition can 
be given as 


Ga 11 

k =6.13X10°°exp\ *7 a) 
Thermo kinetics 

The heat rates determined using equation (11) have 


been compared with the experimentally observed 
heat rates and the results are shown in Figure 5. A 


In k*, min’ 
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Figure 4 Pseudo rate constant plot for the adia- 
batic thermal decomposition of flash composition 
(KNO;. S: Al; 53 : 17 : 30). 
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close examination of Figure 5 indicates that there 
is good agreement between the experimental and 
predicted values. It can be observed now that the 
kinetics obtained in this study are highly reliable. 


Process safety 


Pyrotechnic mixtures are vulnerable to thermal 
hazards. ARC data are used for determining the 
ceiling temperature for processing, handling and 
transportation of hazardous materials. Accordingly 
the practice adopted is that the process/handling 
temperature should be 100°C below the onset 
temperature observed in ARC.”* This rule has been 
in practice in the process chemical industry for the 
safe and successful operation of process plants, 
storage systems and transportation. On these 
considerations, in the case of flash composition the 
ceiling temperature should never exceed 91 °C. 
Although there is no possibility of reaching this 
temperature during normal mixing’ and packing 
processes of flash composition, this temperature 
can be achieved under situations like heat radiation 
from a neighboring area or ignition from unknown 
sources. During such abnormal situations the flash 
composition mixture is vulnerable to hazard. 


Further, impact and friction sensitivities can also 
lead to triggering of explosive decompositions. 
There is no direct correlation available between 
thermal, impact and frictional sensitiveness, either 
to predict one from the other or to predict which of 
these forces can come together to trigger a thermal 
explosion. We hypothesize that an impact or 
frictional stimulus brings about a thermal stimulus 
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Figure 5 Kinetic modeling for the adiabatic 


thermal decomposition of flash composition 
(KNO;. S: Al; 53 : 17: 30). 
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for the flash composition to undergo thermal 
explosion. Under severe impact or friction stimuli, 
thermal stimuli can occur immediately, and this 
can lead to a catastrophic thermal explosion. 
Irrespective of the nature of the stimulus, explosion 
occurs through a thermal mechanism only. This 
means that, for the current flash composition, an 
impact or any other stimulus can only initiate the 
thermal mechanism by providing the minimum 
threshold energy needed/necessary to raise the 
reaction temperature of 191°C as this has been 
observed experimentally as the onset point for 
thermal explosion in ARC. Therefore it is possible 
to relate the mechanical form of energy to the 
threshold energy (AE) observed in the ARC. This 
provides a means of suggesting a predictive 
correlation in such explosive systems. The degree 
of explosivity also depends on the other factors 
such as chemical components, percentages of 
those components, compactness, particle size and 
shape and other environmental conditions. 


Conclusions 


The ARC studies of the flash composition mixture 
confirm that the mixture is vulnerable to thermal 
hazard if exposed above 191 °C. This temperature 
can be achieved under situations like direct heat 
radiation from neighboring areas or ignition 
from unknown sources or through other ignition 
stimuli like impact or friction. For the first time, 
a lower onset temperature of 191 °C for thermal 
decomposition is recorded in this study. It was 
shown that the observed onset temperature can 
be achieved through ignition stimuli like friction 
or impact. Thus this study can offer a better 
explanation for the accident triggering mechanism 
in fireworks factory during the summer months in 
southern India. The kinetic study reveals that there 
is good agreement between the experimental and 
predicted heat rate values. The Arrhenius kinetic 
constants reported in this study are reliable. 
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FIREWORKS — PRINCIPLES AND 
PRACTICE 4" Edition 
Rev R Lancaster MBE 
Review by Christopher Pearce 
Chairman — British Pyrotechnists 
Association 


The name of Ronald Lancaster has been 
synonymous with the British Fireworks Industry 
for the past forty years or so and his firm, Kimbolton 
Fireworks, has been responsible for many of 
the major displays staged in the UK since 1980. 
The fourth edition of Fireworks — Principles and 
Practice is the current, and perhaps final, stage in 
the evolution of a book that was first published in 
1972. 


When the first edition made its appearance, there 
was relatively little in the way of published 
material in respect of fireworks manufacture 
(certainly in the UK). The book was therefore 
seen as a natural successor to classic texts such 
as Weingart’s Pyrotechnics (1947). Over the past 
twenty years or so, there has been an explosion 
(sic) in the number of publications dealing with 
both fireworks and pyrotechnics, at both practical 
and theoretical levels. However, Lancaster’s book 
still has its place and the 4" edition is a further 
refinement on its predecessor, which appeared in 
1998, 


Principles and Practice is presented ina traditional 
format, with twenty-three chapters covering all 
aspects of fireworks manufacture, technology, 
history, and display techniques. In common with 
all three previous editions, there are contributions 
from other well known individuals — particularly 
Tom Smith and Darryl Fleming — whose areas 
of speciality (legislation and display practices 
respectively) give added weight and authority to 
the book. Ronald Hall’s name no longer appears 
in the list of contributors, as Ron Lancaster 
significantly revised the chapter on Rockets in the 
previous printing. 


The fourth edition differs from the third in a few 
important respects. There has been a general 
revision and refinement of the text, a significant 
improvement in the numberand quality of diagrams 
and illustrations, and a fascinating introductory 
chapter entitled Sixty Years’ Love of Fireworks 
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— Diamonds in The Sky. This is the author’s 
personal reflection on a lifetime in fireworks; it 
is a delightful and evocative piece of work. To 
say that it is long overdue is rather crass; after 
all, the right time to pen a worthwhile biography 
is in your twilight years, with all the benefit of 
hindsight (although it is not uncommon for modern 
footballers to record their life history well before 
arriving at their thirtieth birthday!). Lancaster 
takes us back to pre-war days in Yorkshire, where 
his passion for fireworks was kindled, and then 
develops the chapter by recalling his time with 
Pains—Wessex before taking the huge step of 
developing Kimbolton Fireworks as an entity in 
its own right. In forty pages or so, the reader 
is given an insight into a little personal, social 
and industrial history — which is of great value. 
Chapter 1, The History of Fireworks, by Roy 
Butler has not changed significantly since the 
previous edition, but is still an excellent survey of 
the development of fireworks manufacture. It is 
supported by some rather nice images of display 
rigging and factory production relating to Brocks’ 
activities in the 1930s. 


Taking Principles and Practice as a whole — from 
the outset, the author’s philosophy was to produce 
a textbook that adopted a ‘broad-brush’ approach. 
As aresult, all of the key areas involved in general 
firework manufacture are covered in reasonable 
depth, without immersing the reader in minute 
detail. | Take the design and construction of 
cylinder shells, for example. This is covered in 
just a few pages within chapter 16, whereas one 
could devote a whole chapter to do justice to this 
specialised area of shell manufacture. Lancaster 
covers the basics more than adequately, and refers 
the reader to more comprehensive works on the 
subject — a sensible approach that is adopted 
throughout. 


Chapter 8 concerning Mixing and Charging is 
a nice example of how Lancaster’s book has 
evolved since its first printing in 1972. While the 
basic content has remained the same, the chapter 
benefits from improved diagrams (some of which 
are now 3D), a series of new photographs (the 
Pains—Wessex images having been replaced by 
equivalent pictures from the Kimbolton factory), 
and the author’s personal opinions on the current 
state of play. Lancaster does not pull any punches 
when it comes to his assessment of contemporary 
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approaches to health and safety. As an ordained 
Anglican minister, he is able to draw on a source 
of ‘higher authority’ and does so from time to 
time, with amusing results. Bemoaning the 
excesses and complexity of modern health and 
safety management, Lancaster quotes from the 
Bible - Luke XI verse 46 “Yes, you Lawyers and 
Pharisees, you load men with intolerable burdens 
and will not put a single finger to lift the load.” 
Quirky inclusions such as this add a touch of mild 
eccentricity, at the same time enhancing its appeal 
as a work of unique style. 


The key ‘technical’ chapters of the book, as with 
previous editions, cover all fundamental aspects 
of fireworks manufacture and deal with individual 
types; rockets, shells, roman candles, fountains and 
so on. Each chapter is comprehensive enough to 
provide a good basic grounding in the functional 
principles, with a strong practical emphasis 
throughout. The framework of the book has not 
changed significantly, in this respect, compared 
with the 3 edition. Basic firework chemistry 
is also dealt with in a clear and concise fashion 
— providing an excellent grounding in the subject, 
borne out of the author’s own experience. 


The 4" edition of Principles and Practice does 
contain a number of errors — mainly textual and 
typographical, many of which appear in the 
glossary. These do not detract from the overall 
quality of the publication, but need to be addressed 
in any future printings. A couple of important 
‘chemical’ errors have also been identified and the 
author has already produced a list of corrections. 


Ron Lancaster has never needed to strive to 
establish his credentials as an authority in the field 
of fireworks. His research and development work 
has won him considerable respect and Kimbolton’s 
achievements in the display forum speak for 
themselves. Principles and Practice has grown 
alongside Ron Lancaster’s reputation; both author 
and publication have now reached full maturity 
and the 4" edition of this already ‘classic’ text will 
consolidate Lancaster’s position as a key figure 
in the history of fireworks manufacture. This 
latest version should be in the library of anyone 
who has a serious interest in the subject — whether 
that be from a technical or historical perspective, 
or both! 
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2nd Workshop on Pyrotechnic Combustion Mechanisms, 
held July 27th 2005 in Fraunhofer ICT Pfinztal, Germany 


Ernst-Christian Koch? and Rutger Webb? 
“Diehl BGT Defence GmbH & Co. KG, Fischbachstr. 16, D-90552 R6thenbach a d Pegnitz 
>TNO, Defence Safety & Security, Rijswijk, The Netherlands 


For the second time the “Workshop on Pyrotechnic Combustion Mechanisms” was held at Fraunhofer ICT 
in Pfinztal, Germany, preceding the combined 32nd IPS and 36th ICT Seminar in Karlsruhe. The Workshop 
had 41 registered participants from Germany, Poland, UK, USA, Israel, Turkey, Austria, Sweden, Finland 
and the Netherlands. 


This event received positive feedback from both participants and presenters. 
The following presentations have already been published elsewhere. 


. M. Eremets, Polymeric Nitrogen, 
° P. Politzer, Computional Analyses of High-Nitrogen C,H,N Compounds 


° S. Cudzilo, Formation of Carbon Based nanostructures by Combustion of Reductant-Halocarbon 
Mixtures 


° D. Ladouceur, An Overview of the Known Chemical Kinetics and Transport Effects Relevant to Mg/ 
PTFE Combustion 


Journal of Pyrotechnics, Issue 23, Summer 2006 Page 69 


The Combustion Properties of Novel High-Nitrogen 
Energetic Materials 


David E. Chavez, Michael A. Hiskey, My Hang Huynh, Darren L. Naud, Steven F. Son and 
Bryce C. Tappan 


Los Alamos National Laboratory 
High Explosives Science and Technology 
DX-2 MS C 920 
Los Alamos, NM 87545 USA 
naud@lanl.gov and hiskey@lanl.gov 


Abstract: High-nitrogen energetic materials based on the tetrazine and tetrazole ring systems have shown 
unique and unpredictable combustion behavior. Unlike traditional energetic compounds, such as 2,4,6- 
trinitrotoluene (TNT) and hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX), which derive their energy by the 
oxidation of the carbon and hydrogen skeletal atoms by the oxygen carrying nitro group, high-nitrogen 
materials typically have large positive heats of formation as their source of energy. This difference in the 
energy source may partly explain why the combustion chemistries of some high-nitrogen materials are 
unusual. 


Using the precursor 3,6-bis-(3,5-dimethylpyrazol-1-yl)-s-tetrazine (BDT), several useful energetic 
compounds based on the s-tetrazine system have been synthesized and studied. A number of these tetrazine- 
based materials have shown to exhibit burn rates with low sensitivity to pressure, namely 3,6-bis(1H-1,2,3,4- 
tetrazol-5-ylamino)-s-tetrazine (BTATz), 3,6-bis-nitroguanyl-1,2,4,5-tetrazine (NO»Tz), the corresponding 
bis-triaminoguanidinium salt (TAG,NO>Tz) and the N-oxides of 3,3'-azobis(6-amino-1,2,4,5-tetrazine) 
(DAATO3.5). A fifth compound of high nitrogen make-up, triaminoguanidinium azotetrazolate (TAGzT), is 
not prepared from BDT, but it also burns at exceptional rates with low pressure sensitivity. 


The tetrazole-based materials, bis-(1(2)H-tetrazol-5-yl)-amine (BTA) and 5,5'-bis-1H-tetrazole (BT), are 
useful high-nitrogen energetic ligands for the preparation of metal complexes. While BTA, BT and their 
salts have been previously shown as possible energetic fuels for low-smoke pyrotechnic applications, some 
recent combustion experiments with the metal complexes of BT and BTA have proved to be even more 
noteworthy. These metal ion complexes have sufficient internal energy that they can burn under an inert 
atmosphere to produce the free metal, usually in the form of high-surface area foams or nano-sized particles. 
This highly unusual, reductive combustion chemistry may lead to efficient and controlled production of 
metal nanofoams. 


The heat of formation (AH, of 3,6-diazido-1,2,4,5-tetrazine (DiAT), a highly energetic and sensitive 
energetic material (most notably to friction, spark and impact), was calculated to be approximately +1100 
kJ mol, or +92 kJ mol-atom™', using an additive method. Depending on the heating rate, DiAT can 
undergo pyrolytic decomposition to produce either carbon nanospheres or carbon nitride nanopolygons. 
With slow heating, leaf-like or rope-like forms of carbon nitride were the predominant products. With faster 
heating, carbon spheres with diameters on the order of 10 to 100 nm were produced. Such nanomaterials 
are of interest to the scientific community for a wide number of industrial applications. 


Keywords: high-nitrogen, tetrazole, tetrazine, combustion, nanomaterials, propellant, foam 


Introduction materials that exhibit unusual combustion 
properties with a variety of potential commercial 
and military applications. These materials have 
two characteristic attributes. First, they are 


At Los Alamos National Laboratory we have 
developed a number of energetic high-nitrogen 
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high in nitrogen and, consequently, have large 
positive heats of formation, which is attributed 
to the nitrogen—nitrogen bonds that are prevalent 
in these systems. ! Second, because the skeletal 
framework contains relatively less carbon and 
hydrogen when compared to traditional explosives, 
the oxygen balance can be more easily achieved. 
Furthermore, these materials possess relatively 
high crystal densities, a critical parameter for the 
enhancement of explosive performance, and are 
often insensitive towards destructive stimuli such 
as impact, friction, and electrostatic discharge. Of 
particular interest to our group is the synthesis of 
high-nitrogen materials based on the tetrazine and 
the tetrazole systems. In this paper we describe 
the combustion properties of a number of these 
materials in three applications, specifically their 
potential use as high-performance propellant 
ingredients, precursors for high-surface area and 
low-density metal foams, and carbon and carbon 
nitride nanomaterials. 


High-Performance Propellant 
Ingredients 


Solid propellants provide a means of converting 
chemical potential energy into useful kinetic 
energy. Although the chemical ingredients of 
propellants are varied and complex, two basic 
components—fuel and oxidizer—burn to produce 
heat and gas, which can expand in a gun barrel 
to push a shell, or flow at supersonic speed in 
a De Laval nozzle to provide thrust. High- 
nitrogen compounds may be the key to meeting 
the advanced performance objectives of next- 
generation propellants. Classical solid propellants 
have not changed fundamentally for nearly three 
decades and can be separated into two general 
classes. Monopropellants contain ingredients that 
have fuel and oxidizer in the same molecule, such 
as the explosives octahydro-1,3,5,7-tetranitro- 
1,3,5,7-tetrazocine (HMX) or nitroglycerine 
(NG). Heterogeneous or composite propellants 
are aggregate mixtures of fuel and oxidizer 
bound together by polymeric binder, which is 
a fuel itself. The ammonium perchlorate (AP) 
and hydroxy-terminated polybutadiene (HTPB) 
matrix is a popular composite motor propellant 
system. Ingredients, such as HMX, aluminum or 
AP of different particle sizes, are often included 
in the composite matrix to enhance or tailor its 
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performance. A serious drawback of AP solid 
propellants, however, is the copious amount of 
toxic hydrogen chloride that is generated as a 
combustion product. One solid rocket booster of 
the space shuttle, for example, generates 240,000 
pounds (109 tonnes) of hydrogen chloride gas. 


High-nitrogen energetic materials offer the 
possibility of significant performance improvement 
over current propellant systems in at least two 
ways. In one, inert and invisible nitrogen gas is 
the major combustion product of high-nitrogen 
materials, which is a clear improvement over toxic 
hydrogen chloride produced by AP composite 
propellants. This will also facilitate the formulation 
of low-smoke, “reduced-signature” propellants, a 
useful feature that increases the defensive posture 
of an attacker. In the second, high-nitrogen 
materials typically have large positive heats of 
formation and generate low molecular weight 
gases, attributes that are very desirable for high 
impulse performance. This is best exemplified 
by the simplified reciprocal relationship between 
the exit gas velocity, V, of a De Laval nozzle and 
the mean average of the molecular weights of the 
combustion gas products, M,,.. 


V «[1/M,.]? 


ave 


Decreasing the overall molecular weight of the exit 
gases by substituting hydrogen chloride, carbon 
monoxide and carbon dioxide, which are typical 
reaction products of HMX/AP composite motors, 
with lighter gases of high-nitrogen compounds, 
such as H, and N,, will increase the exit gas 
velocity, and in turn the impulse of the motor. 


In this section, we compare the burn-rate data 
of five high-nitrogen energetic materials as 
viable candidates for use in high-performance 
propellants. They are: (1) 3,6-bis(1H-1,2,3,4- 
tetrazol-5-ylamino)-s-tetrazine (BTATz), (2) the 
mixed N-oxides of 3,3’-azobis(6-amino-1,2,4,5- 
tetrazine) (DAATO3.5), (3) triaminoguanidinium 
azotetrazolate (TAGzT), (4) 3,6-bis-nitroguanyl- 
1,2,4,5-tetrazine(NQ,Tz), and(5)itscorresponding 
bis-triaminoguanidinium salt (TAG,NQ,Tz). The 
structures of these materials are given in Figure 1. 
Except for TAGZT, all of the above compounds 
are derived from the valuable precursor, 3,6- 
bis(3,5-dimethylpyrazol-1-yl)-1,2,4,5-tetrazine, 
which is prepared in two steps from commercially 
available triaminoguanidine —_ hydrochloride, 
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2,4-pentanedione, and nitrogen dioxide.* The 
synthesis procedures of these five materials have 
been described elsewhere.* 


In Figure 2, the measured burn rate vs. pressure 
profiles are graphically represented for the pressure 
ranges studied for each material.* The linear 
relationship between burn rate (”,) and pressure 
(p) is best described by the empirical equation 
’»= cp", where c is the empirical constant and n is 
the pressure exponent. A large pressure exponent 
(approaching 1 or greater) is often typical with 
traditional high explosive materials and indicates 
that 2nd order gas-phase reactions dominate the 
combustion process. A lower pressure exponent, 
however, is attributed to early rate-controlling 
reactions in the condensed phase, and thus results 
in a burn rate that is insensitive to changes in 
pressure. Low pressure sensitivity (i.e, a low 


pressure exponent) offers advantages in the design 
of gun and rocket propellants, particularly in motor 
stability. Conversely, energetic materials with 
high exponent values are typically avoided for use 
in propellant applications. In all five materials, 
the pressure exponents are reasonably, if not 
exceptionally low (see Table 1). For comparison, 
the pressure exponent for HMX is 0.75 and its 


burn rate at 1000 psi is 1.1 cms !.° 


Nanoporous metal foams 


The synthesis of low-density, nanoporous materials 
has been an active area of study in chemistry and 
materials science dating back to the first synthesis 
of aerogels.° These materials, however, are mostly 
limited to silica, metal oxides (e.g., Al,O3) and 
organic aerogels (e.g., resorcinol/formaldehyde), 
with the only elemental material being carbon, 
arising from the pyrolysis of organic aerogels. 


Ak \\ + NHNH, 
_ 2 
My 
HN aye 
ok H,NHN“ “NHNH, NN 
N7 “Nn dO Oa 2) 
I | aes _/ N—-N 
NL 2N N=N N=N N =( 
= / = 
bi Ny N NX PN Ne UN 
N. NH N-N 0c —N “Nn * NHNH, 
# \ 
x eO Ob 
N—NH atb+ct+d+e = approx. 3.5 H,NHN NHNH, 
BTATz DAATO3.5 TAGzT 
NNO, + NHNH, NNO, 
lk NH, H,NHN NHNH, Jk NH, 
Yoh ion 
N N ne 
HN 5 a HN = a Tae, 
NNO, NQ>Tz NNO, H,NHN~ ~NHNH, TAG.NQoTz 


Figure 1 Five high-nitrogen energetic materials of interest as high-performance propellant additives: 
(a) 3,6-bis(1H-1, 2,3,4-tetrazol-5-ylamino)-s-tetrazine (BTATz); (b) the mixed N-oxides of 3,3 '-azobis(6- 
amino-1,2,4,5-tetrazine) (DAATO3.5); (c) triaminoguanidium azotetrazolate (TAGzT); (d) 3,6-bis-nitro- 
guanyl-1,2,4,5-tetrazine (NO>Tz); and (e) its corresponding bis-triaminoguanidinium salt (TAG,NQ2Tz). 
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Figure 2 Burn rate versus pressure profile for five high-nitrogen energetic materials. Some of the 
materials were pressed with binders (BTATz with 3% Kel-F-800, DAATO3.5 with 5% PVA and 1% TEG), 
while the rest were studied as neat materials. 


Table 1 Sensitivity and burn rate data for five high nitrogen materials. All data are for pure material un- 


less noted otherwise. 


Pressure Exponent, n 
Empirical Constant, c 


Burn Rate at 1000 psi + 


Hs, (cm) Type 12 
DSC Onset (°C) 
Friction (kg) BAM 
Spark (J) 

Density (g cm’) 
AH; (kJ mol ') 

AH; (kJ mol-atom ') 


BTATz 


0.49" 
0.581 
4.6 

32° 
264 
> 36 
< 0.364 
1.76 
+883 
+40.1 


DAATO3.5 TAGzT 
0.275° 0.67 
1.69 0.287 
5.4 4.9 

50 25 
177 195 
2-14 10 
< 0.36 <0.31° 
1.88 1.60 
+690 +1075 
+29.4 +24.4 


NQ,Tz 
0.163 
0.899 
2.0 

65 
228 
> 36 
> 0.36 
1.76 
+389 
+15.0 


TAG,NQ,Tz 
0.366 
0.494 
2.3 

114 
166 
> 36 
> 0.36 
1.61 
+1255 
+22.4 


*BTATz formulated with 3% Kel-F-800 binder. "DAATO3.5 formulated with 5% PVA and 1% TEG. ‘Positive 
impact results ranged widely between 32 and 200 cm so a conservative value of 32 cm was assigned. “11 out 
of 13 samples initiated when subjected to 0.36 J of energy. “Threshold Initiation Limit (TIL) value, which 
represents a probability that an event will happen 3.4% of the time a given stimulus is applied. + 68.0 atm or 


6.89 MPa. 
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Figure 3 Synthesis scheme of Fe-BTA, formed by the reaction of iron(II) perchlorate with 3 equivalents 
of diammonium salt of bis-(1(2)H-tetrazol-5-yl)-amine (BTA). 


A method to prepare ultra-low-density transition- 
metal foams has been developed using a novel 
pyrotechnic technique. The process is simple in 
concept; a transition-metal complex containing 
an energetic, high-nitrogen ligand is allowed to 
combust under an inert atmosphere, such as argon 
or nitrogen. As the complex burns, the metal 


100nm 
ol P 


*. 


cation undergoes chemical reduction to the free 
metal. Because the intrinsic energy of the selected 
complex is high, the combustion is self-sustaining. 
By this method, metal monolithic foams with 
remarkably low densities (0.011 g cm °) and high 
surface areas (258 m’ g ') have been formed. The 
ability to form monolithic metallic nanocellular 


Figure 4 (a) Photograph, 4 mm scale, of iron foam next to unburned pellet of the Fe-BTA complex. (b) 
Scanning electron micrograph (SEM), 10 um scale, of low pressure iron foam showing pore structure of 
roughly I um. (c) SEM, 100 nm scale, of high pressure iron foam showing pore sub-structure of roughly 
20-100 nm. (d) SEM, I um scale, of heat-treated iron foam. 
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Figure 5 SEM, scale 200 nm, of Co foam after heat treatment, foam walls made up of particles and rods 


of cobalt metal. 


porous materials is presently not possible using 
conventional processes, particularly at the densities 
we have observed or with such ease of production. 
To date, we have produced iron, cobalt, copper 
and silver foams from their corresponding high- 
nitrogen ligand complexes. 


The lowest-density metal foams found in the 
literature range from 0.04-0.08 g cm > and were 
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made from magnesium and aluminum. These 
foams, however, contain cells on the millimetre 
length-scale, and have a relatively low surface 
area.’ Recently, a silver sponge was reported 
with pore sizes on the order of a few microns and 
a surface area of 0.5 m* g ', however, no density 
was given.® The synthesis of dendritic platinum 
with three-dimensional foam-like nanostructures 


nm, 
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Figure 6 SEM, scale 100 nm, of Co foam showing nanoporous substructure made up of rice-like particles 


that form foam walls. 
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Figure 7 SEM, 200 nm scale, of heat-treated copper foam. 


was also recently described, but the structures are 
limited to small clusters ranging in size from 6 to 
200 nm.’ 


The first and best characterized metal foam in our 
initial studies was obtained from the ammonium 
tris(bi(tetrazolato)amine)ferrate(IID complex 
(Fe-BTA). This complex is easily formed by the 
reaction of 3 : 1 molar ratio of the diammonium 
salt of bis-(1(2)H-tetrazol-5-yl)-amine (BTA) 
and iron(III) perchlorate (Figure 3). BTA is a 


high-nitrogen material easily synthesized from 
inexpensive sodium dicyanamide and sodium 
azide under controlled addition of hydrochloric 
acid.'° When ignited as a loose powder in air, 
Fe-BTA rapidly combusts with the production 
of orange sparks, which is attributed to nascent, 
elemental iron burning. The combustion of 
cylindrical pellets of Fe-BTA in a combustion 
chamber of inert atmosphere at varying pressures 
results in the formation of the metal foam monolith 


Figure 8 (a) SEM, scale 1 um, of Ag foam. (b) scale 100 nm, closer view of (a) showing surface nano- 


porosity of around 20 nm. 
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Figure 9 Synthesis scheme for DiAT from 3,6-dihydrazino-1,2,4,5-tetrazine (DHT). 


(Figure 4a). Following combustion at 20.4 atm, 
a bimodal pore size distribution was observed 
within the micron and 20-200 nm ranges. 
Interestingly, combustion at higher pressure (72.3 
atm) resulted in a metallic foam containing only 
the 20-200 nm pore size distribution (Figures 
4b-c). Brunauer-Emmett—Teller (BET) surface 
area analysis for the metal foam produced at 20.4 
atm yielded an extraordinarily high surface area 
of 258 mg |. For comparison, a high surface area 
titania aerogel has BET values that range from 
100 to 200 mg '.'! BET surface areas for foams 
produced at higher pressures (ca. 70 atm) range 
from 12 to 17 mg '. Evidently, the high-nitrogen 
ligand acted as a blowing agent on a molecular level 
as Fe-BTA decomposed, liberating decomposition 
gases. Elemental analysis (standard combustion 
technique) and energy dispersive spectra (EDS) 
demonstrated that the as-produced metal foams 
contain up to 50% carbon and nitrogen which 
are mostly removed by heating to 1073 K in an 
inert atmosphere (Figure 4d). While no density 
or surface area measurements were made on the 
heat-treated foam due to small sample size, it is 
apparent from the SEM image that no sintering 
took place (Figures 4d and 5). Thermogravimetric 
analysis (TGA) and EDS indicate that the heat- 
treated foams range from 10% remaining carbon 


Figure 10 SEM images of carbon nanospheres at 
magnifications of 25,000 x (left) and 150,000 x 
(right). 
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(iron foam) to essentially pure (copper foam). 


Metal foam production using this unique 
combustion method is possible with other transition 
metal complexes. Cobalt, silver and copper metal 
foams have been produced from the corresponding 
BTA complexes; however, the morphology of the 
resultant foam is strongly dependent on the metal 
employed. Electron micrographs of the cobalt 
foam showed the two basic pore morphologies 
observed in the iron foam. A third morphology 
consisting of small cobalt grains (ca. 100 nm) that 
aggregate to form the foam walls, which should 
dramatically increase the surface area of the metal 
foam (Figure 6). The copper foam was of a higher 
density and had much more regular, yet larger pore 
sizes of approximately | micron, and after heating 
to remove impurities, displayed interesting crystal 
lattice lines along the foam walls (Figure 7). The 
silver foam was the most dense, and in fact had 
partially fused to form shiny silver beads. Selected 
SEM images of the silver foam that had not fused 
show a foam structure, indicating that monoliths 
similar to those of iron foam could perhaps be 
achieved by optimizing the combustion method 
(Figures 8a—b). This pyrotechnic method shows 
promise as a flexible and simple approach to the 
formation of a wide range of new nanoporous 
metals that are not currently accessible by state-of- 


Figure 11 SEM images of carbon nanopoly- 
gons at magnifications of 25,000 x (left) and 
150,000 x (right). 
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Figure 12 SEM images of carbon nitride, 
C3N, with leaf-like (left) and rope-like (right) 
morphologies. 


the-art nanoscience. Possible applications abound 
in the fields of catalysis, fuel cells, hydrogen 
storage, unique insulation, and electromagnetic 
absorption materials. 


Pyrolysis of 3,6-diazido-s-tetrazine (DiAT) 


The compound, 3,6-diazido-1,2,4,5-tetrazine 
(DiAT), is easily prepared by the diazotization 
of 3,6-dihydrazino-1,2,4,5-tetrazine | (DHT) 
(Figure 9). Typical of organic azides, DiAT is 
extremely sensitive to spark, friction and impact 
stimuli and must be handled with extreme care 
and in small quantities, preferably stored as 
a solution.’* Using the additive method, the 
heat of formation of DiAT was estimated to be 
relatively high, approximately +1100 kJ mol ', or 
+92 kJ mol-atom '.'* The large positive value is 
attributed to its high nitrogen makeup of 85.4% 
by weight. Because of these attributes, DiAT 
was examined as a precursor for the production 
of carbon and carbon nitride nanomaterials by 
pyrolytic methods. 


Elemental carbon is used in many applications, 
including high-density and high-strength carbon 
artifacts, super-active carbon beads of high surface 
area, lithium storage, lithium battery anodes, 


Figure 13 SEM images of higher-nitrogen 
content carbon nitride, C3;Ns, with sheet-like (left) 
and rope and ball (right) morphologies. 
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spherical packing materials for HPLC, hydrogen 
storage applications, and catalysis.'* Because their 
applications significantly depend on the shape and 
size of the particles, much attention is focused on 
methods of preparation and characterization. 


Besides carbon, carbon nitrides are also of interest 
to the materials community due to their novel 
mechanical, optical, and tribological properties, 
including low density, extreme hardness, surface 
roughness, wear resistance, chemical inertness, 
and biocompatibility. These mechanically hard 
materials promise a variety of technological and 
biological applications. For example, they are 
used as biocompatible coatings on biomedical 
implants, battery electrodes, catalytic supports, 
gas separation systems, electronic materials, and 
humidity and gas sensors.'° Applications of carbon 
nitrides are not only governed by the texture and 
size of the particles but also by the relative nitrogen 
content. As a consequence, extensive effort is 
focused on the discovery of precursors along with 
the appropriate methods to increase the nitrogen 
content in carbon nitrides. 


Carbon nanoparticles were prepared by the 
pyrolytic decomposition of DiAT under 
atmospheric air. The reactions were carried out 
in a 50 mL heavy-walled pressure bomb. For 
ease and safety of handling fairly large quantities 
of DiAT, a method of loading the pressure bomb 
was developed. A chloroform solution containing 
0.2 g of DiAT was poured into the bottom half of 
the bomb and the solvent allowed to evaporate 
in the hood. The bomb was assembled and then 
ramp-heated at a constant rate for 2 hours until 
the temperature reached 150°C. The bomb was 
cooled and disassembled and the resulting fine 
product collected from the walls. Elemental 
and SEM analysis confirmed the production of 
carbon nanospheres with diameters ranging from 
25 to 100 nm (Figure 10). When the reaction 
was repeated, but with the heating accelerated 
such that 150 °C was reached within 1 hour, an 
audible pop within the pressure bomb was heard. 
Elemental analysis indicated that the resulting 
pyrolytic product was again mainly carbon. An 
SEM image revealed carbon nanospheres with 
larger diameters, ranging from 50 to 200 nm in 
diameter (Figure 11). 


Nitrogen-rich carbon nitrides were produced 
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from DiAT when two different heating protocols 
were employed. Using the same pressure bomb 
and under ambient air pressure, DiAT (0.3 g) 
was heated to 100°C over the course of 2 h and 
held at this temperature for an additional 4 h. The 
temperature was then increased to 150°C over 
3h and maintained overnight to yield leaf-like 
carbon nitride CN, (Figure 12, left). When DiAT 
was heated to 150 °C over the course of 5 h and 
then held at that temperature overnight, carbon 
nitride C;N, with a rope and ball morphology was 
obtained (Figure 12, right). The carbon nitride 
products from both reactions were confirmed by 
IR spectroscopy and elemental analysis. Carbon 
nitrides with higher nitrogen content were formed 
when the same two heating protocols as describe 
above were used but conducted under a nitrogen 
atmosphere. As measured by elemental analysis, 
the nitrogen content of the reaction products 
increased from C3N, to C;N;. The SEM images of 
these nitrides show sheet-like (Figure 13, left) and 
rope and ball morphologies (Figure 13, right). 


These results demonstrate a novel method of 
producing carbon nanospheres with diameters 
ranging from 5 to 50 nm and nitrogen-rich carbon 
nitrides of varying morphologies. While DiAT 
may be one possible precursor, other high-nitrogen 
compounds with high positive heats of formation 
may provide additional avenues for the synthesis of 
exotic carbon and carbon nitride nanomaterials. 
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Laser-initiated Reactions Of Energetic/Thermitic Composites 
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Abstract: Researchers are attempting to prepare smaller (nano-scale) metal particles, and nano-scale 
thermitic (metal—metal-oxide) composites. When added to energetic compositions, these nano-materials 
could burn during or close behind the shock front produced by an explosive material. Therefore, investigation 
of their combustion kinetics is important, especially when the investigation technique requires only very 
small quantities of material that is initially prepared. This study uses time-resolved emission spectroscopy 
to measure reaction kinetics and mechanisms of micrograms of material initiated by a laser pulse. Results 
from nano-scale aluminum and the thermite-type compositions of Al+Fe,03, Al+MoO, and Al+B,O; are 


presented here. 


Keywords: /aser, nano-scale, composite, kinetics, time-resolved spectroscopy 


Introduction 


Metals such as aluminum have been added to 
energetic material compositions for a long time. 
Traditionally the Al particles have been in the 
20-100 micrometre range, resulting in Al burning 
after (or behind) the shock front produced by the 
energetic material. More recently, researchers 
have been attempting to prepare smaller (nano- 
scale) metal particles, and nano-scale thermitic 
(metal—metal-oxide) composites. These nano- 
materials could burn during or close behind the 
shock front produced by an explosive material.'~ 
* Jt is important to develop a simple method 
to measure reaction kinetics using the initially 
prepared small quantities of these new materials. 
We use time-resolved emission spectroscopy to 
measure reaction kinetics of microgram quantities 
of sample. In this paper we present time-resolved 
reaction data from nano-scale aluminum and 
the nano-aluminum containing thermite-type 
compositions of Al+Fe,O;, Al+MoO;, and 
AI+B,03, and discuss the kinetic rates of these 
reactive materials. 


Experimental Method 


Sample holders were prepared by exposing 5 um 
gold foil to focused 532 nm laser light from a 5 ns 
pulsed ND:YAG laser. The resulting holes were 
on the order of 80-120 um in diameter. Powdered 
samples were pressed onto the gold foil hole using 


Journal of Pyrotechnics, Issue 23, Summer 2006 


a pestle. Average sample thickness was estimated 
at 15 um. 


Mounted samples were returned to the focus 
position of the 532 nm laser and exposed to a 
single 5 ns pulse. Emitted light was collected and 
sent by a fiber optic cable first into a spectrometer 
to disperse the light spatially, then into a streak 
camera to temporally separate the spectra, and 
finally into a CCD to measure the intensity of 
the light at each time and wavelength. The 
spectrometer was centered at 490 nm and has a 
spatial range of ~100 nm. The streak camera was 
set fora 47 us window. The zero for the time axis 
was determined by the image of the laser pulse 
on the streak camera output. Laser pulse energy 
was estimated from a portion of the beam that was 
diverted into an energy meter using a single glass 
slide. For these experiments the average energy 
per pulse was ~220 mJ. 


Experimental results 
Nano Al 


The energetic/thermitic composites in this study 
all use nanometre size Al as their fuel. Therefore, 
it is important to understand what response should 
be expected from the laser initiation of nano Al by 
itself. 120.5 nm Al was obtained from Technanogy. 
This Al had a 4.9 nm thick oxide coating (74.3% 
active Al content). 
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During all the laser initiation experiments we 
performed, an intense, broad background appears 
in the first couple of microseconds due to the laser 
interaction with surrounding air. This background 
reduces after the first couple of microseconds. In 
the case of nano Al initiation, evidence for AlO 
emission can be seen above the background noise 
as early as the first microsecond. By the third 
microsecond the background has reduced and the 
dominant features in the spectrum are due to AlO 
emission. These results are illustrated in Figure 1. 
Note that due to edge effects on the detector, peaks 
below about 465 nm in this wavelength window 
are significantly less intense relative to the rest 
of the spectrum. The AIO signal that appears by 
3 us persists until 37 us. Beyond that time, any 
peak intensity has reduced to the level of the 
background noise. In Figure | each line plot was 
obtained by averaging over spectra within the time 
range shown in the legend. The stick plot refers 
to emission lines obtained from spectrographic 
references.*° The region from approximately 
520-540 nm is reduced in intensity due to a 532 
notch filter placed in the beam path to reduce the 
intensity of the exciting laser, which defines time 
t=0. 


Al + Fe,O, 


A representative set of spectra for an Al + Fe,0; 


Intensity (arb. units) 


440 460 480 500 


Wavelength (nm) 


sample prepared by Professor Edward Dreizin 
(New Jersey Institute of Technology) is shown 
in Figure 2.’ The first peaks attributed to the 
Al + Fe,0; sample appear above the background 
even during the first or second microsecond. A 
full set of peaks is clearly visible beginning from 
4 us, and the average of the spectra from 4 to 
8 us is shown. The majority of peaks appearing 
in this spectrum can be attributed to elemental 
iron. These results illustrate that the iron oxide 
is breaking apart very early in the reaction. Some 
iron lines persist through 30 us, but any evidence 
of them beyond 38 us is lost in the background 
noise. 


By 30 us there is some evidence for the appearance 
of AlO peaks. The maximum intensity for these 
peaks occurs during 30 to 34 us. A time average 
over this interval is shown in Figure 2. As can be 
seen, the AlO peaks are not well-formed. This 
may be due in part to overlap with the iron lines 
that are still in the spectrum. 


An Al + Fe,O; sample was also obtained from 
Professor Michael Zachariah’s group (University 
of Maryland).* This sample also exhibited Fe 
peaks within the first couple of microseconds after 
exposure. However, by the 9 us peaks from AlO 
emissions become visible and are more intense 
than the Fe peaks. The strongest Fe peaks persist 


520 540 


Figure 1 Time resolved emission spectra from the laser initiation of 120.5 nm Al. The line plots were 
obtained by averaging the spectra within the time ranges shown in the legend. The stick plot was 
generated from emission lines for AlO obtained from a spectrographic reference. 
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440 460 480 500 


Wavelength (nm) 


Al+Fe,O, 38-42us 
Al+Fe,O, 30-34us 
— Al+Fe,0, 20-24us 


— Al+Fe,0, 4-8ys 
? AIO 
* Fe 


Figure 2 Representative spectra for reaction of Dreizin’s Al + Fe,O3. Each Al + Fe,O; line spectrum is 


an average over the time range shown. 


for at least 20 us. The AIO peaks are visible until 
35 us. Beyond 35 us there are no spectral features 
distinguishable above the background noise. Time- 
averaged spectra from the reaction of Zachariah’s 
Al + Fe,O; sample can be seen in Figure 3. 


Al + MoO, 


Figure 4 displays a representative set of spectra 
for the reaction of Al + MoO; (also prepared by 
Dreizin’s group). During the first 8 us of the 


Intensity (arb. units) 


440 460 480 500 


Wavelength (nm) 


reaction there is evidence of elemental Mo above 
the background. This is illustrated by the average 
of the spectra from 2 to 7 us shown in Figure 4. 


By 9 us, weak AIO peaks begin to emerge. These 
peaks gain in intensity, reaching a maximum 
around 22 us. The time average of the spectra 
from 20 to 24 us is shown in Figure 4. By 40 us 
there is still some evidence for AlO peaks, but 
their intensities are barely above the background 
noise. 


Al + Fe,O, 30-34ys 
—— Al+ Fe,0, 20-24us 
— Al+ Fe,0, 9-13pys 
— Al+ Fe,0, 4-6ys 
? Alo 
* Fe 


Figure 3 Representative spectra for reaction of Zachariah ’s Al + Fe,O3. Each Al + Fe,O; line spectrum 


is an average over the time range shown. 
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Wavelength (nm) 


Al+MoO, 40-44us 
— Al+MoO, 20-24ys 
—— AI+Mo0O, 2-7us 
? Alo 
T Mo 


Figure 4 Representative spectra for reaction of Dreizin’s Al + MoO3. Each Al + MoO; line spectrum is 


an average over the time range shown. 


Al+B,0, 


Al + B,O3 samples were prepared at the South 
Dakota School of Mines by Jan Puszynski’s 
group.” The laser initiation of Al + B,O, produces 
a more complicated set of spectra with time. 
For the previous samples, the dominant spectral 
features visible in this wavelength regime could 
be attributed to one species or two species that 
were reasonably well separated in either time 
or wavelength. In the case of Al+B,O; it is 


possible that three species are emitting within the 
same time window. As can be seen in Figure 5, 
during the first 8 us there is evidence for AlO 
emission, but there are other peaks unaccounted 
for by AlO. In fact, by the 9-13 ps range, the AlO 
peaks become less intense relative to the rest of 
the spectrum. The most plausible candidates for 
the extra spectral features seem to be the oxides 
of boron. During the 4-8 us time range evidence 
for BO is apparent. There is also evidence for 


440 460 480 500 


Wavelength (nm) 


Figure 5 Results of the laser initiation of Al + B,O; 
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Table 1 A/O appearance and persistence data from four of the samples studied. 


Sample 

120.5 nm Al 1 
Al+Fe,O; (Zachariah) + 
Al+Fe,O; (Dreizin) 25 
Al+MoO,; (Dreizin) 7 


peaks of BO,. The BO, peaks appear to increase 
in intensity and the BO peaks decrease in intensity 
with time. This suggests that B,O; is not reduced 
all the way to B, but rather, starts to oxide again 
after the initial formation of BO. Apparently the 
oxidation process for boron is competing with that 
of aluminum. 


Discussion 


Table 1 shows a comparison of the AlO appearance 
and persistence for four of the samples from this 
study. The Al+B,O; sample is not included 
because the competing boron oxide signals make 
it difficult to extract information on the AIO signal 
by itself. It is important to note that all of these 
experiments were performed in air. Emission 
lifetimes would therefore not be limited due to 
a lack of available oxygen. The significance of 
the first column relates to the time required for 
the given material to begin releasing energy from 
the reaction of aluminum with oxygen. A smaller 
time in the first column means that the energy 
will be available earlier. The second column 
corresponds to the end of energy release due to 
the formation of AlO. The Al+ Fe,O; sample of 
Zachariah’s group showed a slightly delayed AlO 
emission time as compared to the bare nano Al 
sample. Most likely there is some time required 
to break apart the constituents, which contributes 
to the delay. However, the Dreizin Al + Fe,O; 
sample shows a significantly different behavior. 
The time to first appearance of AIO is over 20 ps 
longer for Dreizin’s material. The probable cause 
for this change comes from the difference in Fe,0, 
content. It is evident even from visible inspection 
that the Fe,O; content in Dreizin’s sample is much 
higher than that for Zachariah’s. Dreizin’s sample 
has an orange-brown color resulting from excess 
iron oxide, whereas Zachariah’s sample is gray 
(indicative of the aluminum). The higher Fe,O, 
content is also illustrated in the persistence of 
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1 appearance of AIO (us) 


Final appearance of AIO (1s) 


41 
ey) 
38 
41 


neutral Fe emission in the spectra. For Zachariah’s 
sample the Fe emission is significantly weakened 
relative to the AlO emission by 9 us, while for 
Dreizin’s sample some Fe emission persists 
beyond 30 us. Such a strong Fe emission may 
also mask the original appearance of AIO in the 
spectra, allowing for the possibility that the AlO 
emission begins at the same time for both samples 
but the emission intensity for Dreizin’s iron 
signal overwhelms the weaker AIO signal early 
in the experiment. For both Al + Fe,0,; samples 
it is clear that the FeO; breaks apart early in the 
reaction. This suggests that this oxide can provide 
a source of oxygen to the aluminum for reaction. 


The Al + MoO; sample shows an initial appearance 
time for AIO similar to that of Zachariah’s 
Al+Fe,0; sample. The initial background is 
higher for the Al + MoO; sample, which may mask 
any possible earlier detection of AlO. Elemental 
molybdenum appears above the background, 
which verifies the release of oxygen. 


For the Al+B,03 sample, the AlO signal is 
convoluted with signals arising from the oxidation 
of boron. While the original B,O; is most likely 
breaking apart initially, any released boron appears 
to be oxidizing at a rate comparable to that of 
aluminum. The various oxide spectral features 
begin to appear around 4 us and last until about 
30 us. 


Similar experiments have been performed on RDX 
with this same experimental method, except a 
backlight was used to track changes in absorption 
because the reaction products of RDX do not emit 
in the visible range available in these experiments. 
A two-stage process was evident in the RDX 
reactions. The first stage was the formation of 
dark intermediate products (evidenced by an 
increase in absorbance). The second stage was the 
formation of clear gaseous products (evidenced by 
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a decrease in absorbance). The time to maximum 
absorbance for reacting RDX should correspond 
to the end of the primary reaction zone. The time 
to maximum absorption for RDX is typically only 
3 us. The time for the entire reaction of RDX 
to reach completion is about 13 ts. Comparing 
RDX reaction times with those of the samples in 
Table | reveals that only the nano aluminum by 
itself shows evidence of AlO formation within the 
time frame of the primary reaction zone of RDX. 
Zachariah’s Al + Fe,O; is close to being within 
this zone at 4 us, followed by Dreizin’s Al + MoO, 
(7 us). However all of the samples in Table 1, 
with the exception of Dreizin’s Al+Fe,O0;, would 
be able to contribute energy from AIO formation 
to the late time (gaseous) reaction of RDX. 


Conclusion 


This study reports on the kinetic and chemical 
processes occurring after the laser initiation in 
air of 120.5 nm Al, Al + Fe,O0;, Al + MoO,, and 
Al+B,03. All of the samples display spectral 
features resulting from the combustion of 
aluminum. AIO emission appears earliest and 
persists the longest for the nano aluminum sample. 
The Al + Fe,O; and Al + MoO; samples both show 
emission from the neutral metals of the original 
oxide (Fe and Mo, respectively) first, followed by 
the appearance of AlO. This suggests that the metal 
oxide is breaking down, which provides oxygen 
for resulting aluminum combustion. However, in 
these experiments excess oxygen was available 
from the surrounding air as well. The Al + B,O, 
showed some evidence for AlO emission, but 
the AlO spectral features were overlapping and 
competing with those of BO and BO). 
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Abstract: In most cases pyrotechnic mixtures are constituted from granular components. A theoretical 
study of such a granular system describes the temperature and concentration evolution in the energetic 
material by using hot spots as source terms for temperature and particles. The progress of the reaction 
is mainly influenced by particle properties which are size, melting and evaporation, gasification and 
surface reactions. In a first step the influence of the particle size ratio between fuel and oxidizer particles 


is investigated. 
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Introduction 


The work presented in this paper is part of some 
ongoing efforts to model combustion and ignition 
processes of propellants and pyrotechnics.” ’ The 
aim of this work is not a perfect simulation of a 
certain process but to investigate single effects and 
their consequences on a system and the method 
applied is to make parameter studies resulting in a 
better understanding of the principal combustion 
mechanism. For a parameter study it is necessary 
to perform a greater number of calculations. There- 
fore models and codes are needed which are fast 
enough to do this in a reasonable time. As a conse- 
quence the models have a reduced complexity but 
high performance. This paper concentrates on a 
model, which is called Hot-Spot Model. 


Model Description 


From aphysical point of view combustion processes 
are mainly dealing with heat transfer: it means heat 
is generated and consumed by different processes 
at different locations. Pyrotechnics are usually 
solid materials. Therefore the model description 
starts with the partial differential equation of heat 
transfer in a solid, where convective effects can be 
neglected (equation 1, below). 


OT[x,t] a _, Of%,t] 
= AT = 
a ee [x7] oo 


There are several analytical and numerical methods 
to solve the differential equation. Here Green’s 
method' is chosen because, if the appropriate 
Green’s function for the homogeneous problem is 
known, it only has to be integrated with the source 
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term of the differential (equation 2, below). 


_ | 
T[x]=JG, (exer) SE aa 
Pp 


The Green’s function for the above differential 
equation in three dimensions is a Gaussian-like 
function (equation 3, below). 


The great advantage of this method is that 
numerical integration is a much faster and a much 
more stable process than differentiation. 


The next step is to choose the source terms. Two 
types of heat sources may be considered. One is 
the combustion enthalpy of the chemical reaction 
of the material. The simplest reaction is of Oth 
order described by an Arrhenius term. 


E 
Dreae = Treo @ 


The other one is the ignition energy which is 
introduced in the system at a certain moment to 
start a reaction. This is the hot spot and it can be 
described by a Gaussian function or a Dirac delta 
function (equation 5, overleaf). 


If there are only hot spots, integration can be 
done analytically. But this is not possible with 
the Arrhenius term for the chemical reaction. 
Therefore numerical integration is needed. The 
algorithm used for this purpose mainly consists of 
three steps. The first step is to generate an initial 
temperature profile resulting from the initially 
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given hot spots (equation 6). 


In the second step the progress of the chemical 
reaction for a small time step At is calculated 
(equation 7). The third step is to calculate 
the heat diffusion for the same time step At by 
integrating the convolution of the temperature 


profile and the Green’s function (equation 8). 


Steps 2 and 3 are repeated for every time step. Step 
1 can be included as often as new hot spots occur. 


From the above description the following features 
can be mentioned: the model is a transient one 
and it describes the heat generation by a chemical 


Figure 1 Temperature profile for a burning particle at the right corner. 
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reaction and the heat transport by conduction. 


Fuel and oxidizer particles as sources of matter can 
also be described as hot spots. Including diffusion 
they can react as soon as gaseous fuel and oxidizer 
come into contact. The model does not include 
any convection or radiation. Implementing the 
model into computer code, we developed several 
versions of the model for different purposes. One 
version includes phase transition. 


To run the program three types of parameters are 
necessary: 


e Different spatial and temporal distributions 
of hot spots: energy, size and number of hot 
spots 


e Material parameters: density, heat capacity, 
heat conductivity, diffusion coefficient 


° Reaction parameters: maximum tempera- 
ture, Arrhenius-parameters (frequency factor 
and activation energy) 


The calculations result in two-dimensional temp- 
erature and concentration profiles for every time 
step, and in heat output, pressure and chemical 
rate over time. 


Results 
Single particle combustion 


First implementation is a one-dimensional version 
of the model. It describes a single particle in 
a reactive atmosphere. As spherical symmetry 
is assumed one spatial dimension is enough to 
describe evaporation or gasification of the particle 
and succeeding combustion. The phase transition at 
the particle surface is implemented as an artificial 
high concentration. No heat transfer to the interior 
of the particle is included. Although these are very 
strong simplifications there are some interesting 
results. 


Figure 1 gives the temperature profiles in a 3 D 
plot. The front edge is the spatial axis and the left 
side is the time axis. The evaporating particle is 
placed at the lower right corner. As no heat transfer 
into the particle is included its temperature is set 
to zero. The surrounding gaseous oxidizer has a 
temperature of 500 K that is enough to evaporate 
a small portion of the particle and to start the 
reaction. In this Figure the evaporation ends when 
temperature reaches the time axis on the right side. 
The flame is marked by the region of the highest 


0.08; » ¢ 
® t 
0.06 = t, = t, _ t 
—_ Parabolic Fit 
0.04 


0.2 0.3 0.4 


Figure 2 Burning time t, = t, — t,; vs. initial radius and parabolic fit. 
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Flame Position 


Figure 3 Flame position as location of the highest temperature vs. time. 


temperature. 


Calculations with different initial radii leads 
to corresponding evaporation times ¢, for the 
particles. If ignition time ¢; is defined as the time 
to reach the flame temperature then the burning 
time 4, = 4, — ¢, can be fitted by a parabolic function 
(Figure 2). This is well known as the 7°-law from 
theory® as expressed in equation 9, overleaf. 


Figure 4 Series of temperature profiles for 
ignition with hot particles. 
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t, = Prey T 
° 2-4, -In(B+1) 


(9) 


On the other hand if the flame position as the 
location of the highest temperature is taken from 
the same calculations Figure 3 results. As can 
be seen the flame moves with time. There is no 
constant distance to the particle surface and the 
flame exists much longer than the particle. Simple 
expressions to describe these findings like the 7’- 
law can not be found in theory, but measurements 
on burning droplets show a similar behavior.’ 


Ignition with hot particles 


For the investigation of ignition processes with 
hot particles a full three-dimensional version 
was developed. A reactive material is hit by hot 
particles which are used as heat sources to start 
the chemical reaction of the reactive material. 
Figure 4 displays a short series of temperature 
profiles. From left to right it shows the initial 
state when the particles hit the reactive material. 
Then the heat spread and chemical reaction of the 
material produces additional heat. Finally a closed 
reaction front is formed. 


The amount of heat produced and the resulting 
temperature can be used as a measure for the 
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Figure 5 Normalized conversion and conversion 
rate during ignition with hot particles. 


progress of the reaction. In Figure 5 it is called 
normalized conversion. The time derivation of 
this quantity is the conversion rate. The peak at 
the beginning of the curve of the conversion rate 
results from the fact that at the beginning every 
particle has its own reactive sphere. Later these 
spheres overlap and the total reactive area reduces 
(see Figure 4). Then the reaction of the material 
goes to a steady state burning. 


Reacting particles 


In most cases pyrotechnic mixtures are constituted 
from granular components, i.e. fuel and oxidizer 
are reactive particles. In the framework of this 
model they can be described as sources of material 
which then react and produce heat. In addition 
to the temperature concentration profiles are 


Figure 6 Zemperature and concentration profiles for eight fuel and eight oxidiser particles of same size. 
A thermal hot spot at the center ignites the material. From left to right: temperature, fuel concentration, 


oxidiser concentration, product concentration. 
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Figure 7 Conversion for 1, 2 and 8 particles of fuel and oxidizer of same size (from left to right) and dif- 


ferent distances between the particles. 


calculated. Figure 6 gives a short series of profiles 
for temperature, fuel concentration, oxidizer 
concentration and product concen-tration. The 
reaction of the particles is started by a thermal hot 
spot in the center. 


For the results displayed in Figure 7 the number 


of particles and the distance between the particles 
were varied. The size was the same for all particles. 
As expected the best burn out is achieved for the 
smallest distance between the particles. For more 
than one particle there is a step in the curve. 


Usually the particles of the oxidizer and the fuel are 


Figure 8 Temperature and concentration profiles for fuel and oxidizer particles of different size. A ther- 
mal hot spot ignites the particles. From left to right: temperature, oxidizer concentration, fuel concentra- 


tion, product concentration. 
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Figure 9 Amount of oxidizer, fuel and product for different O/F-ratios vs. time (from left to right: 7:2, 


7:4. 7:7) 


of different size. Figure 8 shows an example with 
a series of temperature and concentration profiles 
where large oxidizer particles are surrounded by 
smaller fuel particles. 


Varying the initial amount of fuel leads to different 
oxidizer-fuel ratios. Figure 9 gives the curves of the 
total amount of oxidizer, fuel and product vs. time 
for three different O/F ratios. The best O/F ratio 
depends on the type of the chemical reaction. 


In Figure 10 the heat output of the reaction is given 
dependent on the energy of the initial thermal hot 
spot. As can be seen the influence of this parameter 
is less important if the energy is high enough to 
start the reaction. 


As the examples show the hot-spot model is a 
useful tool to make parameter studies. Parameters 
that can be varied are number, size and distribution 
of particles and some material properties (c,, A, p, 
D). Simple chemistry is included. The computer 
code runs fast and is very stable because there is no 
need to calculate derivations numerically. Further 
steps will be to include phase change also in the 
3 D version and more complex chemistry. After 
validation a lot of calculations have to be done 
studying the influence of the different parameters. 
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